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ABSTRACT 

The ALESS survey has followed-up a sample of 122 sub-millimeter sources in the Extended Chandra Deep 
Field South at 870/im with ALMA, allowing to pinpoint the positions of sub-millimeter galaxies (SMGs) to 
^0.3 arcsec and to find their precise counterparts at different wavelengths. This enabled the first compilation 
of the multi-wavelength spectral energy distributions (SEDs) of a statistically reliable survey of SMGs. In 
this paper, we present a new calibration of the MAGPHYS SED modelling code that is optimized to fit these 
ultraviolet-to-radio SEDs of z > 1 star-forming galaxies using an energy balance technique to connect the 
emission from stellar populations, dust attenuation and dust emission in a physically consistent way. We 
derive statistically and physically robust estimates of the photometric redshifts and physical parameters (such 
as stellar masses, dust attenuation, star formation rates, dust masses) for the ALESS SMGs. We find that 
the ALESS SMGs have a median stellar mass M* = (8.9 ±0.1) x 10*** Mq, median star formation rate SFR 
= 280 ±70 Mq yr“*, median overall V-band dust attenuation Ay = 1.9 ±0.2 mag, median dust mass M^asi = 

(5.6 ± 1.0) X 10* Mq, and median average dust temperature Tdust — 40 K. We find that the average intrinsic 
spectral energy distribution of the ALESS SMGs resembles that of local ultra-luminous infrared galaxies in 
the infrared range, but the stellar emission of our average SMG is brighter and bluer, indicating lower dust 
attenuation, possibly because they are more extended. We explore how the average SEDs vary with different 
parameters (redshift, sub-millimeter flux, dust attenuation and total infrared luminosity), and we provide a 
new set of SMG templates that can be used to interpret other SMG observations. To put the ALESS SMGs 
into context, we compare their stellar masses and star formation rates with those of less actively star-forming 
galaxies at the same redshifts. We And that, at z ~ 2, about half of the SMGs lie above the star-forming main 
sequence (with star formation rates three times larger than normal galaxies of the same stellar mass), while half 
are consistent with being at the high-mass end of the main sequence. At higher redshifts (z — 3.5), the SMGs 
tend to have higher star formation rates and stellar masses, but the fraction of SMGs that lie significantly above 
the main sequence decreases to less than a third. 

Keywords: galaxies: ISM - galaxies: evolution - sub-millimeter: galaxies, ISM. 


1. INTRODUCTION 

The first sensitive sub-millimeter bolometer camera (the 
Sub-millimeter Common User Bolometer Array, SCUBA) 
unveiled a new p o pulation of galaxie s over a decade ag o 
jSmail et alJ 1199^ [Hughes et all 119981 : [Barger et"^ 119981) . 
with very large sub-millimeter fluxes (> 1 mJy at 850/rm), 
which were named sub-millimeter galaxies (SMGsfl. These 
galaxies were later identified to be typically at high red- 
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* In this paper we will refer to our sources as SMGs using this purely 
observational definition i.e. our sources are SMGs because they were selected 
to be bright in the sub-millimeter band. In this context, this definition does not 
carry any further assumptions on the nature and intrinsic physical properties 
of these sources. 


shift, and often very faint or comple tely undetected at 
(rest-frame) optical wa v elengths (e.g. ISm ail et al.l 120021 
Dannerbauer et al.l l2002t [Walter et alJ 120121 : ISimpson et alJ 
2014ft . The large sub-millimeter fluxes of SMGs imply very 
large dust infrared luminosities (Ldu.st > 10*^ Lq), which are 
likely powered by intense star formation (exceeding several 
times 100 Mgyr"*). E ven though the n umber density of 
these galaxies is low (e.g. lWeifl et alJl200^ . such highly star¬ 
forming galaxies are thou ght to be the progenitors of lo¬ 
cal el liptical galaxies (e.g. ISmail et alJlT002t ISimpson et all 
12014ft . and the intense gas consumption in the se objects may 
be lin ked to AGN growth and feedback (e.g. iHopkins et alJ 
12008ft . Therefore, understanding these objects i s crucial to 
trace the evolution of todays massive galaxies (see lBlain et alJ 
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I2002t ICasev et^l2014l for reviews). 

A main limitation in characterizing sub-millimeter galax¬ 
ies in terms of redshift and physical properties such as stellar 
mass, star formation rate and dust attenuation has been iden¬ 
tifying their counterparts at shorter wavelengths, due to the 
large beams of the single-dish sub-millimeter discovery ob¬ 
servations. Modern (sub-)millimeter interferometers such as 
the IRAM Plateau de Bure Interferometer (PdBI) and now the 
Atacama Large Millimeter Array (ALMA) are allowing us to 
pin-point the position of SMGs with unprecedented accuracy, 
which allows, for the first time, to reliably identify their coun¬ 
terparts at optical, infrared and radio wavelengths. 

The single-dish LABOCA 870-/im survey of the Extended 
Chandra Deep Lield South (ECDL-S), LESS, is the largest 
and deepest contiguous survey ever performed at that wave¬ 
length, and identified 12 6 sub-millimeter sources with 870- 
fim fluxes above 4 mJy (IWeiB et al.ll2009h . 122 of the 126 
LESS sources were followed-up at the same wavelength, with 
unprecedentedly high sensitivity and spatial re solution en¬ 
abled by ALMA in the Cy cle 0 program ALESS (iKarim et al.l 
1201 3t Hodge et al.l l2013h . The high resolution enabled by 
the ALMA interferometer allowed for a de-blending of multi¬ 
ple sources that were previously identified as a single source 
due to the large beam of single-dish observations, and to 
pin-point t he location of the detected SMGs to within ^ 
0.3 arcsec dHodge et al.l l2013h . Using these observations, 

IHodge et al.l (12013h demonstrated the limitations in identi¬ 
fying the counterparts of sub-millimeter sources in previous 
single-dish studies. The ALESS observations showed that be¬ 
tween 35 and 50% of the detected LABOCA sources are actu¬ 
ally multiple sources blended in the large beam, and that 45% 
of counterparts were missed by previous statistical methods 
that relied on higher-resolution 24-/im and radio observations. 
These observations make ALESS the first statistically reliable 
survey of SMGs, which allows for a complete and unbiased 
multi-wavelength study of the properties of this galaxy popu¬ 
lation. 

The ECDL-S field contains a wealth of ancillary data from 
X-rays, to optical, near - infrar ed, far-infrared and radio wave¬ 
lengths. ISimpson et al.l (1201 4h measured aperture photometry 
of the ALESS counterparts (based on the ALMA positions) 
in up to 19 bands from the f/-band to the Sp/fzer/IRAC 8-/rm 
band. They used the observed optical spectral energy distribu¬ 
tions (SEDs) to compute photometric redshifts, stellar masses 
and dust attenu a tion o f the sources. At longer wavelengths, 
ISwinbank et al.l (l2014ft estimated the infrared fluxes of the 
ALESS SMGs by using the ALMA positions as priors to de¬ 
blend the lower-resolution Herschel photometry, thus building 
the full infrared SEDs from 24 pm to 500-pm, also includ- 
ing the ALMA flux at 870pm and the radio flux at 1.4 GHz. 
ISwinbank et al.l (120141) fit these SEDs using templates to de¬ 
rive the first estimates of total infrared luminosity, dust tem¬ 
perature and star formation rates of these galaxies. 

The ALESS sample and the studies described above allow 
us to build, for the first time, the complete spectral energy dis¬ 
tributions from (rest-frame) ultraviolet to the radio of a com¬ 
plete, statistically unbiased survey of SMGs with robust coun¬ 
terparts. However, so far most studies have analysed the (op¬ 
tical) stellar and the infrared (dust) SEDs of SMGs indepen¬ 
dently from each other (with the exceptioniMichalowski et al.l 
1201 Ol who used GRASIL (ISilva et al.l 119981) templates to 
model the full SEDs of SMGs). Typically, photometric red- 
shifts, stellar masses and dust attenuations are derived from 
optical data, while dust luminosities, dust masses and star for¬ 


mation rates are constrained from far-infrared/sub-millimeter 
data. Moreover, often the models used to analyze the optical 
SEDs are models that were calibrated mostly using ‘normal’ 
local galaxies which are not likely to be as dust-obscured and 
actively star-forming as high-redshift SMGs. This is problem¬ 
atic because stellar age, dust and p hotometric redshift s uffer 
degeneracies (as pointed out by, e.g. [Dunlop et al.l^2007^ . and 
at this high optical depth regime stellar mass-to -light ratios are 
very difficult to constrain ( as discussed by e.g. lHainline et al.l 
1201 Hilsimpson et ^120141) . 

In this paper, we interpret consistently the full ultraviolet- 
to-radio spectral energy distributions of the ALESS SMGs in 
terms of photometric redshift, stellar content, star formation 
activity and dust p roperties using an up dated version of the 
MAGPHYS model (Ida Cunha et alJl2008l) . The new version 
of the code extends the SED parameter priors to the high- 
redshift, high-optical depth and actively star-forming regime, 
mainly by including new star formation histories, a more gen¬ 
eral prior for the dust attenuation parameters, and the effect of 
intergalactic medium absorption of ultraviolet photons at high 
redshift, allowing us to make a detail exploration of the pa¬ 
rameter space for our sources. Additionally, our new infrared 
emission models include a wider range of possible dust tem¬ 
peratures and a simple radio emission component. The new 
SED-fitting code used in this paper also allows for the photo¬ 
metric redshift to be left as a free parameter when fitting the 
observed SEDs. This allows us to explore how well we can 
constrain the physical properties and photometric redshift of 
SMGs when using all available information from the ultravi¬ 
olet to the radio. This paper is organized as follows. In ^ we 
describe the multi-wavelength data available for the ALESS 
SMGs. In ^ we describe the MAGPHYS model, which com¬ 
putes the stellar and dust emission of galaxies in a consistent 
way, and the modifications to the model that were made to bet¬ 
ter explore the parameter space of high-redshift galaxies. In 
0 we present our method to fit the multi-wavelength SEDs of 
the ALESS sources using MAGPHYS, discuss the photometric 
redshifts obtained using our approach, and discuss the prop¬ 
erties of the most optically-faint ALESS sources. In ^ we 
analyze the physical properties related to the stellar content 
and dust properties of our sources obtained with MAGPHYS. 
We discuss the average intrinsic spectral energy distributions 
of the ALESS SMGs, and how the ALESS SMGs compare 
with ‘main sequence’ galaxies at their redshifts in ^ Our 
conclusions are summarized in ^ 

Throughout this paper, we use a concordance ACDM cos¬ 
mology with Ho = 70 km s“' Mpc“*, Ha = 0.7 and Q„ , = 0.3. 
Unless otherwise stated, we adopt a IChabrieij (120031) initial 
mass function (IME) and AB system magnitudes. 

2. MULTI-WAVELENGTH DATA 

Of the 122 LESS sources that were followed-up at higher 
resolution with ALMA as part of the ALESS program, several 
were resolved into m ultiple counterparts , producing a total of 
131 SMG detections (iHodge et al.ll2013l) . 99 of t hese ALESS 
SMGs are considered to be the most reliable bv iHodge et alJ 
(1201 3l) . based on the fact that they are significant detections 
(with S/N> 3.5) inside the ALMA primary beam of good- 
quality maps - these are the sources in the MAIN ALESS 
catalog. The remaining 32 sources (SUPPLEMENTARY cat¬ 
alog) are found to be less reliable bv iHodge et akl (120131) . and 
are therefore are n ot included in the MAIN ALESS sample 
(iKarim et alJl2013l show that up to 30% of SMGs in the SUP¬ 
PLEMENTARY catalog are likely to be spurious). 
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We choose to include only the 99 SMGs from the MAIN 
ALESS catalog in our analysis, because they are the most 
statistically reliable sample of SMGs and therefore allow for 
a complete and unbiased study of this population of galax¬ 
ies. The positions of these SMGs are known to ^ 0.3 arcsec 
thanks to the high-resolution ALMA interferometric observa¬ 
tions, which allows us to identify the counterparts and mea¬ 
sure the fluxes of the SMGs at other wavelengths, using deep 
ancillary data in the optical, infrared and radio. 

2.1. Optical and near-infrared 

To sample the rest-frame stellar emission of our galaxies, 
we use the (aperture) photometry of the ALES S counterparts 
based on the ALMA positions) compiled by ISimpson et al.l 
20141) from archival ground-based and SpitzerllKAC obser¬ 
vations of the ECDL-S; this provides photometry in up to 
19 bands for each galaxy. Most of the UBVRIzJHK pho - 
tometry comes from the MUSYC survey (iTavlor et al.ll2()09h . 
which is supplemented with (deep er) t/-band data from the 
GOODSWIMOS imaging survey (iNonino et al.l 1200^ . In 
the near-infrared, additional deep J- and fCi-band data from 
the ESO-VLT/HAWK-I survey (Zibetti et al., in prep.) and 
the CLHT/WIRC AM Taiwan ECDLS NIR Survey (TENIS; 
[Hsieh et ani2012l) . The SpitzerflRAC data at 3.6, 4.5, 5.0 and 
8.0 /im comes from the Spitzer/IRAC MUSYC Public Legacy 
in the ECDLS survey (SIMPLE; iDamen et al.ll20111) . Lrom 
these image s, seeing- and aperture -matched photometry was 
measured bv ISimpson et al.l(l2014l) across all Alters by center¬ 
ing 3 arcsec apertures on the ALMA position (with a possible 
shift of < 0.5 arcsec) and applying aperture corrections based 
on the total flux of a composite PSL. The uncertainties are 
computed by measuring the flux in 3-arcsec apertures placed 
in random blank regions of the sky in each image (and assum¬ 
ing that sky noise dominates the fl ux uncertainty); full details 
are given in ISimpson et al.l (l2014l) . 

77 of our ALESS main catalog sources are detected (i.e. 
their flux is 3a above t he background noise ) in at least four 
optical/near-IR bandfl ISimpson et al.ll20f4l apply systematic 
offsets to the measured photo metry to optimize SEP fits to 
the HYPER-Z code templates (iBolzonella et al.ir2000l) . Since 
these offsets are SED-model/fitting-procedure dependent, we 
chose to use the raw photometry measured in ISimpson et al.l 
(l2Q14h without including the offsets (in practice this means 
that we subtract the offsets from the photometry quoted in 
their Table 2). In order to account for possible zero point 
offsets and/or mis-matches between the photometry at differ¬ 
ent wayelengths, we add a magnitude error in quadrature to 
the quoted errors t hat is proportional to the photometric off¬ 
sets determined bylSimpson et al.l (l20I4l) . Lor the photomet¬ 
ric bands where ISimpson et al.l (120141) apply smaller offets, 
we conseryatiyely add an additional 0.10 mag to the error in 
quadrature. 

2.2. Mid-infrared, far-infrared and radio 

We use the mid-infrared, fa r-infrared and rad i o flux es of the 
ALESS sources compiled in ISwinbank et al.l (120141) . which 
w e briefly summa r ize he re. 

ISwinbank et al.l (12014t) exploit the Spitzer/MlPS 24-/rm im¬ 
ages of the ECDLS publicly ayailable from the Ear-Infrared 

^ We refer to these SMGs as the ‘optically-bright sample’ and to the re¬ 
maining 22 SMGs which are detected in fewer than four optical/near-IR 
bands as the ‘optically-faint sample’. 


Deep Extragalactic (LIDEL) suryey, and extracted a cata¬ 
log of about 3600 sources down to a 5cr depth of 56/rJy. 
The Very Large Array (VLA) radio obseryation s at 1.4 GHz 
come from catalogue of the ECDLS described in iMiller et al.l 
JMl . The ECDL-S was imaged with Herschel at 70, 100 
and 160p.m with P ACS as pyt of t he PACS Eyolu t ionary 
Probe (PEP) suryey (iLutz et al.l201 ll) . ISwinbank et al.l (120141) 
obtained the PACS fluxes for the ALESS S MGs by match¬ 
ing th em with the PEP catalog extracted by iMagnelli et al.l 
(120131) . HerschellSBYRE imaging at 250, 350 and 500/im is 
ayailable f rom the Herschel M ultitiered Extragalactic Suryey 
IHerMES: lOliyer et al.ll201^ . Due to the large beam of the 
SPIRE obseryations, it is ch allenging to mea s ure th e photom¬ 
etry of indiyidual sources. ISwinbank et al.l (120141) obtained 
the SPIRE fluxes of the ALESS sources using a de-blending 
method based on positional priors from the ALMA maps to¬ 
gether with the MIPS 24-pm and radio source catalogs. 

We note that, since the resolution of the optical/near- 
infrared (~ 1-2 arcsec) and far-infrared (^ 15-25 arcsec) 
are so different, it is not possible to measure the photom¬ 
etry in an identical manner across all wayelengths. How- 
eyer, in all cases, we extract ‘to tal fluxes’: in the optical/near- 
infrared, [simpsoneral] (120141) achieye this by applying aper¬ 
ture corrections based on the total flux of a composite PSP 
determined usin g point sources in the images, as mentioned 
in Section 12.11 and described in detail in section 2.2.1 of 
ISimpson et al.l (120141) . In the far-infrared, ISwinbank et al.l 
(I20I4I) de-blend the Herschel maps by using a PSP at 
each wayelength for each galaxy in the prior catalog, and 
they demonstrate (using simulations) that the total fluxes of 
sources are recoyered using this method. Thus, we expect the 
combination of far-infrared and optical/near-infrared here to 
be as reliable as possible giyen the ayailable data. 

3. THE SED MODEL 

In this paper we aim to model the full obseryed spectral en¬ 
ergy distributions of the ALESS SMGs from the rest-frame 
ultrayiolet to the radio in a phy sically consistent way . To do 
so, we use the MAGPHYS codJi (Ida Cunha et al.ll200^ . which 
relies on an energy balance technique to consistently com¬ 
bine the emission by stellar populations with the attenuation 
and emission by dust in galaxies. In this section, we summa¬ 
rize the main ingredients of MAGPHYS and present an updated 
yersion, which extends the parameter space of the models in 
order to include also properties that are more likely applicable 
to high-redshift SMGs. Specifically, we extend our star for¬ 
mation history and dust optical depth priors (^and ^3.1.2b . 
and we add intergalactic medium absorption in the ultrayiolet 

(EES. 

The stellar population properties of SMGs haye been chal¬ 
lenging to constrain, mostly because of limited (rest-frame 
optical/near-infrared) obseryations ayailable, and also be¬ 
cause preyious analysis haye not modelled the full SEDs 
in a physically consistent way while sampling properly the 
full parameter space of star formation histories, metallici- 
ties, and dust content. Preyious studies that haye attempted 
to deriye the stellar masses of high-redshift SMGs by mod¬ 
elling their rest-frame ultrayiolet to near-infrared SEDs haye 
used limited parameterizations of the star form ation history 
(iHainline et al.ll20TlT: iMichalowski et al]l2012^ . These stud¬ 
ies demonstrate that determining the mass-to-light ratio of 
SMGs is particularly challenging because of the large de- 

^ Publicly available at www.iap.fr/magphys. 
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Figure 1. Distribution of some key physical properties in our stochastic library of models: (a) mass-weighted age (defined as in Eq.|3); (b) mass-to-light ratio in 
the //-band; (c) specific star formation rate; (d) overall V-band attenuation (defined as in Eq.|7); (e) luminosity-averaged dust temperature (defined as in Eq.[^; 
and (f) redshift. We note that the shape of distributions of all properties shown here (except redshift) results from the priors that we chose for other parameters 
of the model that are used to compute them and therefore it is natural that they are not flat (see more details on the chosen priors and how the parameters are 
computed in jOt . 


generacy between stellar age and dust reddening, which is 
exacerbated in the very dust-obscured regime; also, in this 
regime, a large fraction of the mass may be completely dust- 
obscured in optically-thick regions. One of the main goals of 
our SED modelling is to properly take into account these de¬ 
generacies and uncertainties in order to get robust likelihood 
distributions for the stellar masses of the ALESS SMGs. To 
do so, we consider a wide range of SEHs, dust attenuations 
and metallicities, and employ a Bayesian approach to con¬ 
strain the physical parameters from observations. In addition, 
by modelling the infrared dust emission consistently with the 
emission by stellar populations in the rest-frame ultraviolet 
to near-infrared, we expect to better constrain the attenua¬ 
tion by dust and hence get more insight into obscured pop- 
ulations that are required to p ower the observed dust emission 
(Ida Cunha et alJ[2^|2010a|) . 

In the next sub-section, we describe the stellar emission 
models used to fit the UV to near-IR emission from our galax¬ 
ies. These models are build using parameter priors that we 
optimized to get the best possible stellar mass, dust attenua¬ 
tion, mean stellar ages and star formation rate estimates and, 
at the same time, understand uncertainties and degeneracies 
that affect these parameters. 

3 . 1 . UV to near-IR: stellar emission 

We build a library of 50,000 stellar emission models that 
is specifically adapted to interpret the emission from high- 
redshift, possibly dusty sources with unknown ages and star 
formation histories. The prior likelihood distributions for 
each parameter described below are intended to extend the pa¬ 
rameter space of the models (compared to the current publicy- 
available MAGPHYS priors which were calibrated on local 
galaxies of ty pically lower infrared lumino sities and star for¬ 
mation rates: Ida Cunha et al.ll2008l l2010bl) . so that they in¬ 
clude higher dust optical depths, higher star formation rates. 


and younger ages (since we expect galaxies to be younger 
at high redshift). However, it is important to note that we 
avoid biasing our estimates by also including a large fraction 
of models extending to low dust optical depths, low star for¬ 
mation rates and older ages, and allowing our Bayesian fitting 
constrain the most likely values. By doing a Bayesian anal¬ 
ysis and building the posterior likelihood distribution of each 
parameter (based on comparison with the observed pho¬ 
tometry), we ensure robust estimates and specially robust con¬ 
fidence ranges that take into account parameter degeneracies 
and uncertainties due to lack of/poor observational constraints 
in c ertain wavelength ran ges (see more detailed discussion in 
e.g. lda Cunha et al]l2008ll . 

We compute the s tellar emission of each model as in 
Ida Cunha et al.l (l2008l) : essentially, the luminosity per unit 
wavelength emerging at at time t from a model galaxy is ex¬ 
pressed as; 

Lr(0= f t/f'4'(f-f')/fV,Z)exp[-fA(f')], (1) 

Jo 

where lf^^{t',Z) is the luminosity emitted per unit wavelength 
per unit mass by a simple stellar population (SSP) of age f' 
and metallicity Z, T'(f-f') is the star formation rate evolu¬ 
tion with time (i.e. the star formation history), and T\(t') is 
the effective absorption optical depth seen by stars of age t'. 
We compute the emission by SSPs using the s pectra l popu- 
l ation synthesis m odels of [Bruzual & Chariot (120031) (using 
a iChabriCT 120031 initial mass function); we adopt a uniform 
prior in metallicity from 0.2 to 2 times solar. We discuss the 
modelling of star formation histories and dust attenuation in 
the next sections. 

3 . 1 . 1 . Star formation histories 
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For each model, we parameterize the star formation his¬ 
tory (SFH) as a continuous delayed exponential function of 
the form: 

i/’cont (0 oc 7 ^f exp(- 7 f), ( 2 ) 


where t is the model age (i.e. the time since the onset of 
star formation) and 7=1 /tsf is the inverse of the star for¬ 
mation time-scale in Gyr“*. This form of the star formation 
history for high-redshift galaxies is motivated by iFee et all 
(I 2 OI 0 I) . There is growing evidence (based on observations 
and on galaxy evolution models) that, for high redshift galax¬ 
ies, the SFHs are likely to have a variety of shapes, but overall 
they should be rising with time instead of declining as the 
‘r-models’ typically used at low redshift (iPacifici et al.ll201^ 
iBehroozi et atllMl l^ha et al.ll20l4 . The SFH parame¬ 
terization we adopt here (Eq. |2]i rises linearly at early ages 
and then declines exponentially with the timescale defined by 
the 7 parameter. For each model, we draw 7 randomly from 
a broad prior distribution (7 can vary between 0.075 and 1.5, 
which corresponds to the peak of the SFH happening between 
0.7 and 13.3 Gyrs after the onset of star formation). We draw 
the age randomly between 0.1 and 10 Gyr. This ensures that 
our model library includes a wide range of possible contin¬ 
uous SFHs that go from essentially r-like models peaking at 
very early ages, to approximately constant star formation rates 
with time, to SFHs that are linearly increasing with time. De¬ 
spite the fact that we cover a wide range of star formation evo¬ 
lution with this prior, the actual form of SFHs is likely to be 
more complex than any analytic al parameterization (see e.g. 
IConrovl2013HSimha et al.l2014ft . To account for stochasticity 
on the star formation histories, we superimpose star formation 
bursts of random duration and amplitude to the continuous 
model defined in Eq. |2] We set the probability of a burst of 
star formation occurring at any random time in the previous 
2 Gyrs to 75%; each burst can last between 30 and 300 Myr, 
and the total mass of stars formed during the burst can have 
any random value between 0.1 and 100 times the mass formed 
by the underlying continuous SFH (Eq.|2]i. This ensures that 
we account for as wide a range of SFHs as possible using 
our simple parameterization (including both starburst-like and 
more quiescently star-forming), which is crucial to sample all 
possible stellar ages and mass-to-light ratios in our analysis, 
and thus get the most robust constraints on these parameters 
for the observed galaxies. 

In this context, the ‘age’ of a galaxy model is simply the 
time when the star formation history starts in our analytical 
parameterization and has no real physical meaning. In order 
to have a more reliable measure of the overall age of the stellar 
population in our galaxies, we define the ‘mass-weighted age’ 
of each model as: 


ageM 




(3) 


where 'k(f-f') is the star formation history of each model, 
which is essentially 'ipcontit -1') (Eq. |2l plus random bursts. 
The value of age^ hence depends not only on the model age, 
but also strongly on the shape of the star formation history. In 
Fig- eta), we plot the distribution of age^ in our model library, 
that results from the different random SFHs and model ages 
included. We note that /?-band light-weighted ages (i.e. the 
ages of the stars dominating the rest-frame /?-band light) for 
the models in our library are typically lower than the mass- 
weighted ages, with age^/age^ = 0.82 ±0.30. 


We define the current star formation rate (SFR) of each 
model as the average of the star formation history T* over the 
last 10 Myr. 


3 . 1 . 2 . Dust attenuation 

We use the two-component model of ICharlot & Falll (120001) 
to describe the attenuation of stellar emission at ultraviolet, 
optical and near-infrared wavelengths. This model accounts 
for the fact that young stars that are still inside their birth 
clouds are more dust-attenuated than intermediate-age and old 
stars in the diffuse ISM. In practice, the effective optical depth 
in Eq.[T]is described as: 


where 


rx(t') 


-BC , 7 ISM 
-ISM 


for t' < ffic, 
for t' > fee, 


f = (1 - p) fv(A/5500A)“''3 


(4) 

(5) 


is the effective attenuation optical depth of dust in stellar birth 
clouds, and 


“ = p fy ( A/5500A)“‘’'^ (6) 


is the effective attenuation optical depth in the diffuse ISM. 
In practice, the three free parameters of our dust attenua¬ 
tion model are the lifetime of stellar birth clouds t^c, the 
effective (i.e. angle-averaged) V-band optical depth seen by 
stars younger than t^c in the birth clouds fy, and p, the 
fraction of Ty seen by stars older than fee- We change 
the priors of these parameters compared to previous studies 
(Ida Cunha et al.ll?008l 12010b^ in order to reflect our uncer¬ 
tainty of stars/dust geometry and optical depth in high-redshift 
SMGs and to ensure that we cover a broad parameter space. 
We allow for dust attenuation to reach higher values both in 
the stellar birth clouds and the diffuse ISM, to account for 
the fact that S MGs might be heavily optically-thick (such as 
local UFIRGs; Ida Cunha et akll^lOal) . The prior probability 
distribution for fy is constant between 0 and 6 , and declines 
exponentially to a maximum of 20 ; for p, we set the prior 
probability distribution to be a Gaussian centered at 0.25 with 
0.10 standard deviation. We let the lifetime of birth clouds 
fee vary between 5 and 50 M yr - typically this was fixed at 
10 Myr (iCharlot & Fallll20()0h : since this is an unknown pa¬ 
rameter and the ISM conditions in SMGs might be different 
than in moderately star-forming local galaxies (e.g. higher 
densities, stronger stellar winds from the starburst), we vary 
this unknown parameter so that the spectral fitting solution 
marginalizes over a wide range of possible values. 

Even though our dust attenuation prescription is described 
by the three parameters described above, for simplicity and 
to allow for a more direct comparison with other studies, we 
also define, for each model, the resulting overall F-band dust 
attenuation as r em 

AF = -2.51og-^, (7) 

where Ly'’® is the emitted (i.e. observed) V-band luminosity, 
and Ly“‘ is the intrinsic (i.e. dust-free) V-band luminosity of 
each model. The priors of fy, p and fsc described above result 
in the prior distribution for Ay plotted in Fig.lTfd). 


3 . 1 . 3 . IGM absorption 

We include absorption by the intergalactic medium (IGM) 
in our model SEDs, which strongly affects the rest-frame ul¬ 
traviolet emission from high-redshift galaxies. To do so, we 
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use the absorption prescription from iMadaul (Il995h . which in¬ 
cludes Lyman series line blanketing and Lyman-continuum 
absorption. To account for different opacities along different 
lines of sight (which reflect different distributions and prop¬ 
erties of absorbers along the line of sight), we draw the IGM 
effective absorption optical depth of each model from 

a Gaussian distribution center ed at the mean IGM effective 
absorption absorption given in iMadaul (Il995h . ^JGM.Madau^^^^ 
with a standard deviation of 0.5. 

3.2. Dust emission and radio component 

We model the emission by dust from rest-frame mid- 
infrared to millimeter wa velengths using f our m ain dust com¬ 
ponents as described in Ida Cunha et al.l (120081) : (i) a poly¬ 
cyclic aromatic hydrocarbons (PAHs) empirical template; (ii) 
mid-infrared continuum from hot dust; (iii) warm dust in ther¬ 
mal equilibrium; and (iv) cold dust in thermal equilibrium. 
The dust components in thermal equilibrium emit as modified 
black bodies [u^B^iT)], with emissivity index /3 fixed at 1.5 
for the warm components and 2.0 for the cold components. 
The equilibrium of warm dust in stellar birth clouds is 
uniformly distributed in the range 30 - 80 K, and the equi¬ 
librium temperature of cold dust in the ‘diffuse ISM’ 
is uniformly distributed in the range 20-40 K. The relative 
contributions of these dust emission components to the total 
infrared luminosity are also free parameters of the model. 

We note that the detailed PAH emission in SMGs may be 
different from our adopted template, and there may even be 
variations of the PAH emission between SMGs, since the 
strength of different features depends on factors s uch as the in¬ 
tensity of th e radiation held and metallicity (e.g. lDraine & Lil 
1200111200^ . Previous obs ervations of the mid-infrared emis- 
sion from z ^ 2 SMGs by iMenendez-Delmestre et al.l (120091) 
show that their average mid-infrared spectrum is similar to 
that of local starbursts, which im plies that our PAH template 
is at le ast a good approximation. IMenendez-Delmestre et al.l 
(|2009|) also And some variation of the ratios of different PAH 
features from SMG to SMG that we are not able to reproduce 
using our fixed template, however due to the lack of data sam¬ 
pling the (rest-frame) mid-infrared emission for our SMGs in 
detail, we would not be able to constrain that part of the spec¬ 
trum, even if we included variation of the PAH emission in our 
modelling. Furthermore, while the shape of the PAH emission 
in our model is Axed, the contribution of PAHs to the total dust 
luminosity is allowed to vary, therefore the models include a 
range of overall PAH emission strength. While this parameter 
itself can be hard to constrain in the absence of photomet¬ 
ric data sampling the rest-frame mid-infrared, thanks to the 
MIPS 24-/rm data (combined with the Herschel and ALMA 
data which constrain the total infrared SED), we are able to 
loosely constrain the PAH contribution for our galaxies. This 
contribution is typically less than 10% of the total infrared lu¬ 
minosity, and therefore uncertainties related to the PAH tem¬ 
plate should not affect the energy balance significantly. 

Even though our infrared SEDs contain multiple temper¬ 
ature components, for simplicity, and to make comparisons 
between galaxies easier, we define an average, luminosity- 
weighted dust temperature for each model as: 


dust components, respectively, and is the fraction of to¬ 
tal dust luminosity contributed by the diffuse ISM compo¬ 
nent (which are all free parameters of the model). and 
are the fractional contribution and temperature of warm 
dust to the cool ISM component, and are fixed at = 0.0 7 
and = 45 K (see model details in Ida Cunha et al.ll2008l) . 
By varying the temperatures and relative contributions of the 
warm and cold dust to the total dust emission, we obtain the 
distribution of average dust temperatures in our model library 
plotted in Fig. [He)- 

To fit the observed radio fluxes of the ALESS SMGs, we 
also add radio emis sion to our SED s by fo llowing the simple 
method described in ida Cunha et al.l (12013l) . which essentially 
uses the radio/far-infrared correlation and fixed slopes for the 
ther mal and non-thermal radio emission for each model (see 
also [Dale & Heloull2002^ . To summarize, we compute the 
radio emission as the sum of a thermal (free-free emission) 
component with spectral shape L* oc ' and a non-thermal 
component with shape L"* cx We Ax the contribution 

of the thermal component to the 20-cm radio continuum at 
10%, and scale the total radio emission of each model to the 
infrared luminosity using the local radio/far-infrared correla¬ 
tion q i.e. the ratio between the far-infrared and the 1.4 GHz 
flux density. We draw the value of q randomly for each model 
from a Gau ssian prior distr ibution centered at g' = 2.34 (the 
local value: r^n et alJl2001l) . with a cr^ = 0.25, in order to ac- 
coun t for possible scatt er in the radio/far-infrared correlation 
(e.g. Ilvison et al.l 120101) . A lack of strong redshift evolution 
of the radio /far-infrared correlation is supported by several 
studie s (e.g. Ilbar et al.ll2008l: Ilvison et akllioiOt [Urgent et af] 
l2010t [Huang et al.ll2014t iThomson et al. 2014h . and our fits 
(e.g. Fig. 12) show that this assumption is sufficient to repro¬ 
duce the observed 1.4 GHz flux of most ALESS sources de¬ 
tected in the radio. 


AGN contamination. — We note that the SED models used in 
this study assume that the dominant source of dust heating 
is star formation, and do not include AGN emission. This 
should not have a great impact in our analysis of the full 
sample for two main reasons. First, the fracti on of AGN 
in our SMG sample is likely low, as shown by IWang et aH 
dMl, who And that only 10 of the 91 ALESS SMGs in the 
area covered by Chandra are X-ray sources, and 8 of those 
(ALESSOll.l, ALESS017.1, ALESS057.1, ALESS066.1, 
ALESS070.1, ALESS073.1, ALESS084.1, ALESS 114.2) are 
identified as AGN. Second, most of the parameters recov¬ 
ered by MAGPHYS are robust to AGN contamination even 
when the AGN contributes up to 75% of the ultraviolet- 
to-infrared emission, as shown in an independent study by 
[Hayward & Smitfil (l2015l) . The parameter that suffers most 
uncertainty in the case of strong AGN contamination is the 
stellar mass, which may be overestimated by MAGPHYS by 
up to 0.3 dex when the contribution of A GN-heated dust to 
the n ear-/mid-infrared is the highest (see [Hayward & Smitlil 
120151 for more details). 

3.3. Combined ultraviolet-to-radio model SEDs 


Tdust = 


er X T^^ + iT X + X X /, 


( 8 ) 


where and are the fraction of total dust luminosity 
contributed by the (birth cloud) warm and (diffuse ISM) cold 


An important feature of our model is that we combine the 
stellar emission (attenuated by dust) with the dust emission in 
a self-consistent way using a simple energy balance argument: 
that the energy absorbed by dust in stellar birth clouds and 
the diffuse ISM (T®y^, and respectively) is re-radiated 
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in the infrarecfl. Different combinations of physical parame¬ 
ters (star formation histories, dust attenuation parameters) can 
lead to the same absorbed energies and in a model 
galaxy, and these energies can be distributed in the infrared 
range in different ways depending on the parameters regu¬ 
lating the contribution and temperatures of the different dust 
components. In practic e, w e associate each model in our stel¬ 
lar population libr ary (33.11 1 with all the models in the infrared 
spectral library ( ^3.21 l that have similar fraction of the total 
dust luminosity contributed by the diffuse ISM component, 
U = lZ/ILII + L^] (within an eiTor interval Sf^ = 0.15 
which accounts for uncertainties in dust modelling due to ge¬ 
ometry, for example). 

4. MULTI-WAVELENGTH SED EITS 
4.1. Fitting method 

We use the Bayesian method described in Ida Cunha et al.l 
(l2008h to compare our library of SED models described in 
the previous section with the observed photometry in the ul¬ 
traviolet, optical, infrared and radio for each ALESS source. 

The standard version of the MAGPHYS code fits the SEDs at 
fixed redshift. Since spectroscopic redshifts are not available 
for the whole sample of ALESS SMGs (A. Danielson et al., in 
prep.), and the uncertaintie s associated with the optical/near- 
IR photometric redshifts of ISimpson et al.l (12014ll can be quite 
large (when few photometric data points are available to con¬ 
strain the fits), here we leave the redshift as a free parameter in 
our analysis, and effectively test the use of MAGPHYS, for the 
first time, as a ‘photometric redshift code’. This has two ad¬ 
vantages that we discuss in detail later: (i) we can incorporate 
the far-infrared and radio data as additional constraints on the 
photometric redshift (which is particularly important for the 
most optically-faint SMGs), and (ii) we can include the un¬ 
certainties on the photometric redshift in the error bars of all 
other physical parameters in a self-consistent way, since the 
likelihood distributions of the physical properties and redshift 
are computed simultaneously. 

Eor each model in our library, we compute the predicted 
flux in our 28 bands (from the t/-band to 1.4 GHz) for a set of 
100 redshifts drawn randomly from the prior distribution plot¬ 
ted in Eig.dtf), which includes a broad peak at z ~ 2.5 (note 
that if the model age is older than the age of the Universe 
at a given redshift, than that model-i-redshift combination is 
not considered). We use this prior in order to avoid oversam¬ 
pling unlikely regions of the redshift space, however we have 
checked that using a flat redshift prior does not change our re¬ 
sults. At each redshift, the predicted model flux in each band 
is computed by first applying the IGM absorption presc ription 
to the stellar emission at that redshift as described in ^3.1.31 
and then convolving the total (stellar-i-dust emission) model 
SED in the observed-frame with the filter response functions. 
Then, we compare the observed fluxes of our galaxies in all 
the observed bands with the predicted model fluxes by com¬ 
puting the goodness-of-fit for each model in our library 
(upper limits are included by setting the flux to zero and 
adopting the upper limit value as the flux uncertainty). We 

^ As in Ida Cunha et all I2008I) . dust self-absorption is not included. Dust 
might be optically-thick to its own radiation on some lines of sight, which 
is not included in the modelling, since this would require a more complete 
radiative transfer calculation. This should not significantly affect our results, 
because dust self-absorption affects mostly the rest-frame mid-infrared spec¬ 
trum, which is not sampled in detail by our data, and the angle-averaged op¬ 
tical depth values are relatively small. Moreover, dust self-absorption should 
not affect our energy balance to more than a few percent. 


then build the likelihood distribution of each parameter in our 
model (including the redshift) by marginalizing the probabil¬ 
ity of each model P oc exp(-x^/2) over all other model pa ram- 
eters (more details can be found in Ida Cunha et al.l[2008ll . We 
take our estimates of each parameter to be the median of its 
likelihood distribution, and the confidence range as the 16th 
to 84th percentile range. 

As an example, in Fig|2l we plot our best-fit model (i.e. the 
model that minimizes x^) SED obtained with this method, as 
well as the fit residuals in all bands, for the first source in our 
catalog, ALESSOOl.l (the fit results for all the other sources 
can be found in the online version of the journal). For com¬ 
parison, we also plot the SED of the prototypical low-redshift 
ULIRG, Arp220, normalized at the observed ALMA 870- 
/rm flux, at the photometric redshift of ISimpson et al.l ([2013), 
obtained from fitting so lely the I/-band to IRA C fluxes us¬ 
ing the HYPER-Z code (iBolzonella et al.l [2000h . Assuming 
that the SMGs have similar SEDs to that of a typical local 
ULIRG, then if the HYPER-Z photometric redshift is correct, 
the Arp220 SED should approximately follow the observed 
photometry in all the other bands. If this is not the case (e.g. 
ALESS014.1, ALESS015.1 and others) then this means that 
either the photometric redshift is not accurate, or that the SED 
of these SMGs is different from that of Arp220 (or both). This 
is why it is important to fit the whole SED using a wide range 
of possible SED shapes (which we obtain by varying the phys¬ 
ical parameters in our spectral library), and leaving the red¬ 
shift free, as we do in this study. 

In the bottom panels of Fig.|2l we plot the normalized prob¬ 
ability distributions for several parameters (redshift, stellar 
mass, star formation rate, mass-weighted age, average dust 
attenuation, //-band mass-to-light ratio, total dust luminos¬ 
ity, average dust temperature and dust mass), obtained when 
fitting the whole SED, and also when fitting only the ultra¬ 
violet, optical and near-infrared fluxes, and when fitting only 
the mid-infrared to radio fluxes. These probability distribu¬ 
tions show how tightly our method constrains the redshifts 
and physical properties of the ALESS SMGs, and also high¬ 
light the effect of fitting the whole SED simultaneously versus 
fitting only the stellar emission or only the inf rared emission 
as done in previous studies o f these SMGs (e.g. lSimpson et all 
12014t[Citibank et al.ll20l4) . In panels (f), (g), (h) and (j), we 
also plot the normalized prior distribution of each parameter 
in our model library (same as shown in Fig. [T]), in order to 
compare directly between the prior and posterior likelihood 
distribution. The stellar mass, star formation rate, dust lumi¬ 
nosity and dust mass depend on the model scaling, and there¬ 
fore they do not have fixed prior distributions. The stellar 
mass prior is essentially flat, and the other parameters scale 
with stellar mass in a way that is determined by the model 
prior; for example, the SFR will be determined by the stellar 
mass and by the specific star formation rate prior, as plotted 
in Fig. He). 

Our Bayesian fitting method implies that, if a parameter 
is unconstrained (because e.g., of a poor sampling of the 
SED, degeneracies, or large observational uncertainties) the 
posterior likelihood distribution of a given physical parame¬ 
ter should be the same as the prior distribution. When data 
(i.e. multi-wavelength broad-band fluxes) are available to 
constrain a given parameter, the posterior likelihood distribu¬ 
tion becomes different, as each model in the prior distribution 
is now weighted by our well it fits the data, and not all models 
in the library fit the data equally well. For the majority of our 
SMGs, the posterior likelihood distributions of the parame- 
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Figure 2. Example of the MAGPHYS outputs for the galaxy ALESSOOl.l (all other SMG fit figures ai‘e available on the online version of the journal). 
panel: Spectral energy distribution fit of the galaxy. The circles show the measured photometry and the arrows indicate 3f7 upper limits, as described in |2] 
The black line shows our best-fit SED model to all the bands from the ultraviolet to the radio. In the bottom we plot the fit residuals in each band, computed 
as W e also show, for comparison, the SED of the local ULIRG Arp220 scaled to the measured ALMA HlO-pm flux at the redshift found 

bv ISimpson et al.l j2Q14l) by fitting the ultraviolet/optical/near-IR data with HYPER-Z. In the top-left corner, we indicate our median-likeli hood redshift obtained 
by with MAGPHYS (with the 16th-84th percentile range indicated in brackets), and also the HYPER-Z photometric redshift obtained bv ISimpson et alj <20141) 
(with their Icr range in brackets). Bottom panels: (a) joint likelihood distributions of the V-band dust attenuation and redshift; (b) joint likelihood distribution 
of the dust temperature and redshift; the remaining panels show the (normalized) marginalized probability distribution for the redshift (c), stellar mass (d), star 
formation rate (e), mass-weighted age (f), overall F-band dust attenuation (g), ff-band mass-to-light ratio (h), dust luminosity (i), dust average temperature (j) and 
dust mass (k). In all these panels, the black lines indicate the probability distributions obtained when fitting the whole ultraviolet-to-radio SED. The green lines 
and the orange lines show the likelihood distributions we obtain when fitting only the ultraviolet-to-near-infrared (i.e. stellar emission) and mid-infrared-to-radio 
(i.e. dust emission), respectively. In panels (c), (0, (g), (h) and (j) we also plot, in grey, the prior distribution of each parameter, for comparison. 
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Figure 3. Comparison between our median-likelihood redshift estimates and 
the spectroscopic redshifts available for 45 SMGs in the main ALESS sample 
from the zLESS program (Danielson et al., in prep.)- The 19 sources with the 
most robust spectroscopic redshifts are highlighted with squares. The vertical 
error bar shows the 16th - 84th confidence range of our redshift likelihood 
distributions. The points are colour-coded by number of available detections 
in the optical bands (i.e. from U to IRAC- 8 /j.m). 

ters plotted in the bottom panels of Fig. |2] are different from 
the prior distributions, which shows that these parameters are 
constrained by the data (with the width of the distributions 
indicating how well constrained the parameters are, and typ¬ 
ically becoming narrower as more data points are available 
to sample the SED and/or observational uncertainties become 
smaller). We find that the typically least well constrained pa¬ 
rameters in our analysis (which can be seen from the fact that 
the posterior likelihood distributions often resemble the prior 
distributions) are the mass-weighted ages and the dust tem¬ 
peratures. It is not surprising that the ages are not tightly con¬ 
strained given that the rest-frame optical/near-infrared SEDs 
are often not well sampled by the available data, and that 
the stellar age is a very challenging parameter to constrain 
from broad-band data alone (because of age/dust attenuation 
degeneracies, and the fact that young massive stars outshine 
older stars). Regarding the dust temperatures, they are not 
tightly constrained in cases where there are no available Her- 
schel detections sampling the peak of the dust SED (as is the 
case for ALESSOOl.l; Fig.H. 

4.2. Photometric redshifts 

5. COMPARISON OE PHOTOMETRIC AND SPECTROSCOPIC 
REDSHIETS 

In the next section we analyze the properties of the ALESS 
SMGs based on the likelihood distributions of their physi¬ 
cal parameters obtained using the method described above. 
The likelihood distribution of each parameter is obtained by 
marginalizing the probability over all other parameters, in¬ 
cluding the redshift. Here we compare our MAGPHYS-based 
photometric redshift estimates (using the full optical-to-radio 
SEDs) with the spectroscopic redshifts obtained by A. Daniel¬ 
son et a l. (in prep.), and with the photometric redshift esti¬ 
mates in iSimpson et al.l (120141) . in order to investigate the reli¬ 
ability of our redshift probability distributions before moving 
on to the analysis of the intrinsic physical properties of the 
galaxies. 

Spectroscopic redshifts are available for 45 of the main 
sample ALESS galaxies from the zLESS program (A. Daniel¬ 
son et al., in prep.). In Fig. |2 we compare our MAGPHYS 


Figure 4. Comparison between the ISimpson et all 120141) photometric red¬ 
shifts obtained by fitting the U -band to SphzerflRAC photometry using the 
HYPER-Z code IBolzonella et all 1200(1) (zsis), and the MAGPHYS median- 
likelihood redshifts obtained from fitting the full IZ-band to radio SEDs 
(zmagphys)- The error bars show the statistical confidence ranges (for the 
MAGPHYS redshifts, these correspond to the 16th to 84th percentiles of the 
likelihood distribution marginalized over all other model pa rameters). The 
orange symbols correspond to the 22 ‘optically-faint’ SMCJs: ISimnson et all 
120141) do not constrain individual redshifts for these sources, but they esti¬ 
mate a median redshift of 2513 ~ 3.5 for the SMGs detected in two or three 
optical bands, and zsi 3 — 4.5 for the SMGs detected in less than two bands. 
The inset shows an histogram of the relative difference between the two red¬ 
shift estimates for the 77 ‘optically-brigh t’ SMGs for which we have individ¬ 
ual estimates from lSimpson et Hi 120141) (green symbols). 
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Figure 5. Normalized stacked redshift probability distribution of all 99 
SMGs in our sample (in black). The green histogram shows the contribu¬ 
tion by the 77 ‘optically-bright’ sources (with four or more detections in the 
optical/near-infrared range), and the orange histogram shows the contribution 
by the remaining 22 ‘optically-faint’ sources. In grey , we plot, for compar- 
ison, the stacked probability distribution obtained by ISimpson et all 12()14f) 
using the HYPER-Z code. The vertical lines show the medians of the distribu¬ 
tions (also indicated in the top right corner). 

photometric redshift estimates with those spectroscopic red¬ 
shifts. We hnd a generally good agreement between our pho¬ 
tometric redshift estimates and the zLESS redshifts. On av¬ 
erage, we hnd a median relative difference of Az = (Zspec “ 
2magphys)/( 1 +2spec) = “0.005, with a standard deviation of 
0.29. Many of the largest outliers have very large photomet¬ 
ric redshift conhdence ranges that are still consistent with the 
spectroscopic redshift. However, there are nine SMGs with 
Az > 0.3 and conhdence ranges that do not include the spec- 
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troscopic redshifts; ALESSOlO.l (robust Zspec), ALESS029.1, 
ALESS087.1 (robust z^pec), ALESS071.1, ALESS084.1, 
ALESS069.1, ALESS080.1, ALESS055.1, ALESS037.2. In 
the case of ALESSOlO.l, ALESS084.1, ALESS069.1, the 
(well-sampled) observed SED seems to be inconsistent with 
the spectroscopic redshift, particularly the MIPS 24-p,m de¬ 
tection tends to constrain the redshift to be lower. In the 
case of ALESS037.2, redshift seems inconsistent with the 
UV/optical upper limits. Eor ALESS071.1 and ALESS055.1, 
our MAGPHYS fit does not optimally reproduce the IRAC 
near-infrared data and there is a large observed excess, so 
probably our fit is not trustworthy. Assuming that the ob¬ 
served SEDs and the spectroscopic redshift are correct, some 
of these outliers could be cases where the energy balance im¬ 
posed by MAGPHYS does not work properly and so our pho¬ 
tometric redshifts may not be reliable. This could happen if 
the stars/dust geometry in the galaxies is very different than 
what is assumed in the models, for example if most of the 
rest-frame UV/optical light comes from a region that is phys¬ 
ically separated from w here the bulk of the infrared/sub-mm 
emission originates (e.g. [Hodge et al.ll2012h . 

In Eig |4] we compare the median-likelihood estimates of 
the photometric redshift of our sources obtained by fitting 
the full ultraviolet-to-radio SED with MAGPHYS, Zmagphys^ 
with previous estimates using the HYPER-Z photometric red- 
shi ft code (which only fits the stellar emission) obtained 
by ISimpson et al.l (l2014l) . zsi 3 . Overall, we find that the 
agreement between the redshift estimates is good and there 
are no strong systematic offsets. The inset histogram in 
Eig. a shows the distribution of the relative difference be¬ 
tween the redshift estimates, which we quantify as Az = 
(zsi3-Zmagphys)/(1 +ZS13)- The median of this difference is 
only 0.008 (which is much smaller than the typical error bars), 
and the standard deviation is 0.32. The largest outliers in this 
comparison are ALESS087.1, ALESS124.1, ALESS061.1, 
ALESS006 1 (likely gravitatio nally-lensed), ALESS083.4 
(flagged in ISimpson et al.l [20141 as possibly having contam¬ 
inated optical/near-IR photometry), and ALESS065.1 (de¬ 
tected only in three IRAC bands in the optical/near-IR range). 
In Appendix |5] we compare our photometric redshifts with 
spectroscopic redshifts available for 45 of our SMGs (Daniel¬ 
son et al., in prep). 

Our method allows us to investigate in detail the advan¬ 
tage of fitting the full observed t/-band-to-radio SED simul¬ 
taneously in deriving the physical properties of galaxies and 
also photometric redshifts. The bottom panels of Eig. |2] illus¬ 
trate the differences between parameter likelihood distribu¬ 
tions obtained when fitting the full SED and when fitting only 
the optical-to-near-infrared emission, or only the infrared-to- 
radio SED. We find that fitting the full SED presents a unique 
advantage particularly for galaxies for which the number of 
detections in the optical/near-infrared is low. One exam¬ 
ple of this is ALESS002.2, which in the optical/near-infrared 
range is detected only in three IRAC bands. The analysis in 
ISimpson et al.l (|2014|) is unable to constrain the photometric 
redshift for this source. When we fit only these three opti¬ 
cal fluxes using our code, we also obtain a very broad like¬ 
lihood distribution for the redshift (and also the other phys¬ 
ical parameters). The joint redshift-Ai/ probability distribu¬ 
tion shows a clear degeneracy between these two parame¬ 
ters (this is a known degeneracy in photometric redshifts and 
also noted bv iDunlop et al.ll2007^ . When we add the infrared 
data to the fit, since we have four infrared fluxes (from Her- 
schellSPlRE and ALMA) constraining the peak of the SED, 


the redshift is better constrained as well as the other parame¬ 
ters, particularly the dust luminosity and Ay are constrained 
via energy balance between the optical and infrared. This 
is slightly alleviated when more detections sampling the dust 
emission are available from Spitzer/MlPS and HerschelfPACS 
(e.g. ALESS014.1, ALESS015.1). When only one infrared 
flux is available (e.g. ALESS015.3, ALESS023.7), the red¬ 
shift is still very hard to constrain, because with only one far- 
infrared flux we cannot constrain the peak of the SED. 

How well each photometric redshift is constrained depends 
on the number of photometric bands where each galaxy is 
detected, and also on its intrinsic redshift and SED. There¬ 
fore, we characterize the redshift distribution of the whole 
sample in Eig. |5] by stacking the marginalized redshift like¬ 
lihood distributions of the galaxies, which naturally ac¬ 
counts for the different redshift uncertainties affecting dif¬ 
ferent sources. The median of this stacked likelihood dis¬ 
tribution is Zphot = 2.7 ±0.10, with a 16th - 84th percentile 
range from 1.6 to 4.1. We note that th is redshift distr i bution 
is consistent with the one derived in ISimpson et alJ (l2014t) 
for the same sample ( as shown in Ei gJUL,, and also with 
other SMG samples (e.g. [Ch apman et al.|2005]:| Aretxaga et alJ 

iTOn^lWardlow et al.lUOl lUSmolcic et akiTOm 

5.1. The optically-faint ALESS sources 

Before we move onto a detailed discussion of the physi¬ 
cal properties of all 99 SMGs in the main ALESS catalog 
(©, here we discuss the 22 SMGs that have the lowest de¬ 
tection rate in the optical/near-infrared range, the ‘optically- 
faint’ sub-sample. These sources are detected in fewer than 
four bands in the study of the optical/near-infrared em ission 
of the ALESS counterparts of ISimpson et al.l (l2014h . Live 
of these sources (ALESS035.2, ALESS055.2, ALESS069.3, 
ALESS087.3, ALES103.3) are undetected in the t/-band to 
IRAC-8/im photometry, as w ell as in the infrared w ith Spitzer 
or Herschel and in the radio (ISwinbank et al.ll2014^ . 

The redshifts and physical properties of these galaxies are 
particularly challenging to constrain due to the lack of data 
sampling the SEDs. Our SED modelling approach allows us 
to attempt to characterize these faint galaxies by combining 
all the (scarce) available data in the optical/near-infrared and 
far-infrared/submillimeter, and by including the information 
available from the upper limits. In Eigs. |5] and |6] we plot 
the stacked likelihood distributions thus obtained for these 22 
SMGs and compare them to the stacked likelihood distribu¬ 
tions of the 77 ‘optically-bright’ sources (i.e. SMGs that are 
detected in four or more optical/near-infrared bands). We also 
include the medians and 68% ranges of these distributions in 
Table □ 

The likelihood distributions of most properties of the 
optically-faint SMGs are consistent with them having simi¬ 
lar properties to the optically-bright SMGs. The main excep¬ 
tions are the likelihood distributions of the redshift and av¬ 
erage V-band dust attenuation Ay. The median of the red¬ 
shift likelihood distribution for the optically-faint SMGs is 
Zphot = 3.7 ±0.1, significantly higher than the median red¬ 
shift of the optically-bright SMGs (Zphot = 2.3 ±0.1). This 
difference in the average redshift of these so urces is consis¬ 
tent w ith the stacking analysis performed in ISimpson et al.l 
(120141) . who concluded that the optically-faint SMGs de¬ 
tected in none or one optical band have an average redshift 

^ We refer to the MAGPHYS-derived photometric redshifts as Zphot from 
now on. 
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Figure 6. Normalized stacked likelihood distributions of different physical parameters for all 99 ALESS SMGs: (a) stellar mass; (b) star formation rate; (c) 
specihc star formation rate; (d) mass-weighted age; (e) average V-band dust attenuation; (f) mass-to-light ratio in the H band; (g) total dust luminosity; (h) 
luminosity-averaged dust temperature; (i) total dust mass. We also plot the contribution by the 77 ‘optically-bright’ SMGs (detected in at least four optical/near- 
infrared bands; in green) and the 22 ‘optically-faint’ SMGs (detected in fewer than four optical/near-infrared bands; in orange). The vertical lines indicate the 
median of each distribution (listed also in Table[T). The likelihood distributions show that many of the parameters of the optically-faint sample are unconstrained 
or consistent with the optically-bright sample. However, the optically-faint sample shows significantly higher average dust attenuation Av than the optically-bright 
sample (which also causes the inferred H-band mass-to-light ratios to be larger). 


of 4.5, while the SMGs with two or three detections have 
an average redshift of 3.5. The average dust attenuation in¬ 
creases from Av = 1.9 ±0.3 for the optically-bright sources 
to Ay = 2.9 ±0.3 for the faint sources. These results suggest 
that the optically-faint ALESS SMGs may be a subset of the 
high-redshift (z > 3) SMG population that have higher dust 
optical depths. Since our analysis suggests that the total dust 
mass and total dust luminosity of the sources is similar to the 
SMGs in the optically-bright sample, then these higher dust 
optical depths are likely caused not because there is more dust 
in these galaxies, but because the stars/dust geometry is differ¬ 
ent. Eor example, a higher Ay can be a result of edge-on view¬ 
ing angle, or of a more compact galaxy where dust column is 
higher; this could be tested by comparing high spatial resolu¬ 
tion sub-millimeter continuum imaging of the optically-bright 
and optically-faint sources at similar redshifts. Our results im¬ 
ply that we could be biasing the redshift distribution of SMGs 
low when the redshifts are based on detected optical counter¬ 
parts, because we could be missing this population of very 
obscured, high-redshift sources that is a non-negligible frac¬ 
tion of all the SMGs in the ALESS sample (22%). Obtaining 


redshifts through the detection of molecular lines and/or [CII] 
in the sub-millimeter is the only reliable way to get the real 
redshift distribut ion of these sources (e.g. [Walter et al.ll2012|; 
iWeifi et al.ll2013h . 

In the following sections, we analyze the physical proper¬ 
ties of the complete sample of 99 ALESS SMGs as a whole, 
i.e. including both the optically-faint and optically-bright 
galaxies. 

6 . THE PHYSICAL PROPERTIES OF ALESS SMGS 

In order to analyze the overall properties of our sample con¬ 
strained with MAGPHYS, and at the same time take the un¬ 
certainties associated with these constraints into account, in 
Eig. |2]we stack the all 99 individual likelihood distributions 
of stellar mass, star formation rate, specific star formation 
rate, mass-weighted age, V-band dust attenuation, //-band 
mass-to-light ratio, total dust luminosity, average dust temper¬ 
ature and total dust mass. We note that in our framework, if 
the parameters of the galaxies were poorly constrained, these 
stacked (posterior) likelihood distributions would resemble 
the prior distributions (Eig. [T]i and would change depending 
on the adopted priors. We tested this by htting all SEDs with a 
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Figure 7. Black histograms: normalized stacked likelihood distributions of different physical parameters for the 99 galaxies in our sample, as shown in Fig.|6] The 
blue histograms show the stacked histograms for the sub-sample of galaxies with Zpi,ot < 2.7 (45 sources), and the red histograms correspond to the sub-sample 
of galaxies with Zpi,o, > 2.7 (54 sources). The medians and 16th-84th percentile ranges of these likelihood distributions are listed on Table[T] 


model library with significantly different priors for the param¬ 
eters, specifically including a larger fraction of models with 
older ages, lower specific star formation rates, higher mass- 
to-light ratios, and a flat redshift distribution. We find that 
even when using these different parameter priors the result¬ 
ing stacked posterior distributions and the average properties 
listed do not change significantly, which indicates that our re¬ 
sults are robust, i.e. not biased by the priors. 

In Table[T] we list the medians and the 16th-84th percentile 
ranges of these distributions for each parameter. These ranges 
reflect the range of values for each parameter found in the 
sample, taking into account the uncertainties in each indi¬ 
vidual constraint. Overall, the population of ALESS SMGs 
seems to have fairly uniform properties, with star formation 
rates, dust luminosities, mass-weighted ages, mass-to-light 
ratios and dust masses within < 1 dex of the sample median. 
The range of dust attenuation and temperature values is rela¬ 
tively large (with 16th-84th percentile ranges 0.9 < Ay < 3.1 
and 33 < Tdust/K < 53), but we note that the individual like¬ 
lihood distributions tend to be broader for these parameters. 
To understand if the overall properties vary between SMGs at 
lower and higher redshifts, we also plot, in Eig. |7] the con¬ 
tribution to the stacked likelihood distributions from sources 
with Zphot <2.7 and fro m sou rces with Zphot >2.7 (the median 
redshift of the sample; ^4.21 l. The average properties in these 


two sub-samples are also listed in Table[T] This shows that the 
spread in properties shown in Eig.|2]is also due to some extent 
to the fact that there is a tendency for SMGs at higher red- 
shifts to have higher star formation rates, dust luminosities, 
dust temperatures and dust attenuation values. Eig. |7] also 
shows that our SMGs have a wider range of stellar masses, 
specific star formation rates and mass-weighted ages in the 
low-redshift bin (z < 2.7), with some indication that the dis¬ 
tribution of these parameters may even be bimodal. This may 
be because, at lower redshifts, we may have relatively more 
massive galaxies with older stellar populations that are also 
bright in the sub-millimeter. At higher redshifts, we may not 
be able to detect massive, relatively more ‘quiescent’ galax¬ 
ies because they have not had time to build up to high stellar 
masses and older ages. Additionally, since the SEDs of the 
highest-redshift galaxies are typically less well sampled in the 
optical than those of the lowest-redshift galaxies, it is possi¬ 
ble that in the higher redshift bin any existing structure in the 
distributions of these parameters may be diluted by the large 
uncertainties. 

6.1. Stellar masses & mass-to-light ratios 

To understand how well we can constrain the stellar mass 
of individual SMGs, we analyze their individual stellar mass 
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Table 1 

Average properties of the ALESS SMGs 




Whole Sample 


Optically Bright 

Optically Faint 

Parameter 

all redshifts 

^phot 2.7 

^phot ^2.7 

all redshifts 

all redshifts 


[99 SMGs] 

[45 SMGs] 

[54 SMGs] 

[77 SMGs] 

[22 SMGs] 

log(M*/MQ) 

10.95t« 

10.75)«:« 

io-95!S:4 

ii.05)jj:5 

I0.85t^ 

log(SFR/M 0 yr-') 

2.45^^ 

2.25)“:| 

2.65t^ 

2 

^•^^-0.5 

2.45)|>:^ 

log(sSFR/yr-*) 

-8.55!«:S 

-8.65t? 

-8.45ti 

-8.55)<>;f 

-8.45)j>:^ 

log(ageM/yr) 

8.35ti 

8.45)|>:? 

8.35ti 

8.35)|>i 

8.35;[]i 

Av 

1 0+1.2 
^•^-1.0 


2 4+1.5 

1 ■+!:!! 

2 9+2.2 

log(M*/LH/M0/L0) 

-0 n+0-4 

-0-i8!gi 

-0.03!)>j 

-0.18!<>j 

0-1+04 

( ^dust /^O ) 

12.55^^ 

12 2S+^-^ 
iz.zd_o_4 

12.65![]:| 


12.45)<;j 

TiujK 

43!10 

38fl|> 

43flj> 

38);‘5 

43)1° 

log(Md„s,/M0) 

8.75t^ 

8.75!«:3 

8.65!)Jj 

8.75)|>i 

8-6+o:4 


Note. — Average properties of the ALESS galaxies inferred from stacking the likelihood distributions of the physical parameters of individual 
galaxies. The parameter values listed are the median of the stacked likelihood distribution of each the whole sample of 99 galaxies (first column) and 
each sub-sample (following columns). The range indicated with each median corresponds to the 16th - 84th percentile of the likelihood distribution. 
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Figure 8. Median-likelihood estimates of the //-band mass-to-light ratio of the 99 ALESS sources plotted against the median-likelihood estimates of (a) the 
mass-weighted age, (b) the stellar mass. On the bottom right-hand corners, we plot the median error bars (defined as the 16th-84th percentile range of the 
likelihood distributions). Each galaxy is color-coded according to its derived V-band attenuation. Ay (a), and its //-band absolute magnitude (b), obtained by 
dividing the stellar mass of each galaxy with its inferred mass-to-light ratio. We note that the typical error bars on our stellar masses are smaller than for the 
mass-to-light ratio because we do not use directly to constrain M*. In MAGPHYS, each model SED is nor malized to 1 Mq, and the stellar mass of a 

galaxy is constrained by finding the SED normalization that best fits the data using all available bands (as described in ida Cunha et al.l2QQ8h . 


likelihood distribution^. The stellar mass confidence range of 
each individual galaxy, given by the 16th - 84th percentile of 
the likelihood distribution (marginalized over all other param¬ 
eters, including redshift), has a median of 0.5 dex. This means 
that given the uncertainties caused by degeneracies between 
parameters, unknown star formation histories, sampling of the 
SEDs, and errors in the photometric data, we constrain the 
stellar masses of individual objects typically within a factor 
of 3 (see also iHainline et all 1201 Ih . The error bar on each 
individual stellar mass estimate depends quite strongly on the 
number of detections in the optical/near-infrared range, which 
is expected, since the more detections are available, the eas- 

® In this paper ‘stellar mass’ refers to the current mass of stars in the galaxy, 
i.e. excluding the mass returned to the interstellar medium. 


ier it is to break degeneracies and constrain the parameters. 
Eor the optically-faint galaxies, detected in fewer than four 
optical/near-IR bands, the median confidence range is signifi¬ 
cantly larger at 0.7 dex (or a factor of 5), while for the galaxies 
with the most ideal sampling of the optical SED (with more 
than 10 bands), the median conhdence range is 0.3 dex wide 
(or a factor of 2). 

The median stellar mass o f the ALESS so urces is 8.9 ± 
0.1 X 10'° Mq (assuming a IChabri^ 120031 slightly 

higher than the value found bv ISimpson et al.l (120141) for the 
same sample (8 ± 1 x 10'° Mq for a Salpeter IME, which cor¬ 
responds to about 5 ± 1 X 10® Mq for a Chabrier IME). The 
16th-84th percentile range of the stellar mass stacked like¬ 
lihood distribution is quite broad, from 1.4 x 10^° Mq to 
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~ 5.6 X 10" Mq, which shows that statistically, there is a 
wide range of stellar mass values that are consistent with the 
available data for these galaxies. 

The stellar mass confidence ranges become on average 
^ 0.2 dex broader if only the t/-band to IRAC photometry is 
used to constrain the fit (as shown in the likelihood distribu¬ 
tions plotted in Fig.|2]l. This difference in confidence ranges is 
particularly pronounced for the optically-faint SMGs, where 
the stellar emission-only fits yield a median confidence range 
of 1 dex. This shows that adding the infrared emission con¬ 
straints in our SED fits can help better constrain even the 
stellar masses of these galaxies. This is mainly because the 
infrared information typically reduces the uncertainty in red- 
shift, and it also reduces the uncertainty in dust attenuation 
Ay (via the constraint on total dust luminosity and energy bal¬ 
ance), which affects directly the mass-to-light ratios as shown 
in Fig.|8la). 

We analyze how the stellar mass and mass-to-light ratio es¬ 
timates from our full SED fits depend on other physical pa¬ 
rameters in Eig. in In Eig. EJa), we plot the //-band mass- 
to-light ratio against the mass-weighted age, which is pri¬ 
marily determined by the star formation history. The mass- 
to-light ratio correlates with age^j^ (and star formation his¬ 
tory); the older the ages, the more stellar mass can be ‘hid¬ 
den’ in old, low-mass stars that contr ibute little to the ligh t. 
This has been discus s ed also by e.g. [Hainline et all ( 1201 ll) ; 
ISimpson et al.l (|2014|) ; iMichalowski et al.l (120141) . The large 
range of mass-to-light ratios with different star formation his¬ 
tories illustrates why we must marginalize over a wide range 
of SFHs in our SED fitting to get reliable stellar mass esti¬ 
mates. We find that the typical mass-weighted ages of our 
sample are less than 1 Gyr (the median for the sample is 
~ 230 ± 140 Myr), which indicates that the stellar mass is not 
dominated by a signific antly old stellar population (with ages 
> 1 Gyr) as argued by IMichalowski et all (l2012h . We note 
that, even though the mass-weighted age is a very hard pa¬ 
rameter to constrain, this result does not seem to be driven by 
our modelling assumptions, since our model library includes 
models with age^ >10^ yr, as can be seen in the prior his¬ 
togram plotted in Eig. [Ha). The crucial d i fferenc e between 
our analysis and that of iMichalowski et’aH (120121) is that we 
include a wide range of star formation histories, ages, metal- 
licities, dust attenuations and redshifts, and when we fit the 
whole ultraviolet-to-radio SED we require self-consistency 
between the stellar and dust emi ssion (in terms of energy bal¬ 
ance, as detailed in Section IXTI i. When doing so, it becomes 
clear that these galaxies must have a significant amount of on¬ 
going star formation to power the large infrared luminosities, 
and even when allowing for models to have a significant frac¬ 
tion old stars (our star formation history prior also includes 
models with early bursts/peaks of star formation, which have 
older mass-weighted ages), that is not what is preferred by the 
fitting in a statistical sense, i.e. the posterior likelihood distri¬ 
butions of the age are not peaking at the oldest ages allowed 
by the prior. If indeed the data preferred older ages, we would 
see that the posterior likelihood distributions would peak at 
log(ageM) > 9 and even towards the edge of allowed ages. 
The relation between mass-to-light ratio and age is broad¬ 
ened by dust attenuation even in the H band; Fig. Eta) shows 
that, at fixed age, the higher Ay, the higher the //-band stellar 
mass-to-light ratio. This shows how important it is to prop¬ 
erly model the effects of dust in the stellar SEDs via our en¬ 
ergy balance technique to better constrain Ay and hence the 
mass-to-light ratios and stellar masses. 



log (SFR / Mg yr"') 

Figure 9. Median-likelihood estimates of the star formation rate against the 
total dust infra red luminosity of the ALESS s ources. The bla ck line shows 
the relation of IKennicuttI 0998^ (scaled to a IChabri^ j2003f) IMF) which 
essentially assumes an optically-thick starburst. 


InFig.EIb), we plot the //-band mass-to-light ratio against 
our stellar mass estimates, where each galaxy is color-coded 
by its rest-frame //-band absolute magnitude (computed by 
dividing the stellar mass of each galaxy by its inferred mass- 
to-light ratio). The median Mh inferred from our SED fits 
is Mh = -22.9 ± 1.2. At fixed Mh, the stellar mass corre¬ 
lates tightly with the mass-to-light ratio, and we find that the 
highest mass galaxies have higher M,/Lh, meaning they have 
possibly a larger number of older stars and/or higher dust at¬ 
tenuation. 


6.2. Star formation rates 

We find a median star formation rate for the ALESS 
sources of 280 ± 70 Mg, vr~^ which is consistent with the 
value inferred by ISwinbank et alJ (120141) using the total in¬ 
frared luminosity (300 ±30 Mq yr“^). The star formation 
rates are consistent because, as shown in Fig. |9l the total 
dust luminosity traces the star formation rate extremely well 
(within ^ 0.1 de x), in the same way as predicted using the 
iKennicuti (Il998l) standard conversion between infrared lumi¬ 
nosity and SFR (plotted as a black line). This is expected 
for very dusty, a ctively star-forming sources for which the 
iKennicutl (119981) conversion is best calibrated (as also dis¬ 
cussed by e.g. ^nnicutt et al.l 1200^ Ida Cunha et akll^lOal ; 


iRowlands et al. 20141) . Fig. 1^ also shows that galaxies with 

older mass-weighted ages typically have slightly higher dust 
luminosity per u nit star formation rate (up to ~ 0.15 dex) than 
predicted by the iKennicutll (Il998l) conversion. This is a con¬ 
sequence of additional dust heating by relatively old stellar 
populations (that are more important in galaxies with higher 
age^), which increases the dust luminosity at fixed star for¬ 
mation rate. 

The very tight correlation between the infrared luminosity 
and the star formation rate that we find is not directly imposed 
in MAGPHYS, since the model al lows for significant he ating of 
dust by old stellar populations (Ida Cunha et alJ[^08l) . which 
could imply larger dust lumino sities per unit sta r formation 
rate than those predicted by the IKennicutl (Il998l) conversion 
or, in other words, the fit would need lower SFRs to repro¬ 
duce the observed dust luminosity. While models with non- 
negligible heating of dust by older stellar populations are al¬ 
lowed in the libraries we use to fit the data, the results indicate 
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Figure 10. Luminosity-averaged dust temperature of our SMGs as a func¬ 
tion of their total dust luminosity. Each SMG is color-coded according to its 
measured ALMA 870-/im flux. In the top left-hand corner we plot the me¬ 
dian error bar on these properties, defined as the 16th - 84th percentile range 
of the likelihood distributions. 


that the observed SEDs (including energy balance) constrain 
our sources to be very dusty, actively star-forming galaxies 
where the dust emission is tracing the SER. 

6.3. Dust properties 

On average, the typical ALESS source has a dust luminosity 
of(3.5±0.8)x 10*2 Lq, adust mass of (5.6± 1.0) x IO^Mq, 
and a luminosity-averaged dust temper ature of 43 ± 2 K. 


Kovacs et al.l 120061: IWardlow et al. 

I201C; Mamelli et al.l 

20121 Svmeonidis et al.l 120131 

Swinbank et al.l 120141 


Casey et al 


2014I) . we hnd a correlation between the total 


dust luminosity and the average temperature of our sources 
(pl otted in Eig. rra. The measured dust t e mpera tures 
in IChanman et alJ (120051) and ISwinbank et al.l (120141) are 
typically lower than our luminosity-averaged tempera¬ 
tures, but this could be simply because of the different 
temperature dehnition and/or the way the infrared SEDs 
are h t ted (as discussed in e .g. lCasevlu012l iHavward et alJ 
1201 2t iMagnelli et al.l |2012|) . Our luminosity-temperature 
correlation is broader tha n that observed for t he SMGs of 
IChanman et al.l (l2005h and iKovacs et al.l (l2006l) because our 
sample includes lower sub-millimeter flux sources thanks to 
the high sensitivity of ALMA. At hxed infrared luminosity, 
the deeper we go in sub-millimeter flux 587o^m, the higher 
dust temperatures (and lower dust mass) we probe. This 
means that SMGs only have very cold dust temperatures 
(compared to local galaxies of similar infrared luminosity) if 
we consider the sources detected at th e highest sub-millimeter 
fluxes (see also IMagnelli et aO2012h . 

We hnd a weak correlation between the total dust luminos¬ 
ity (and the average dust temperature) and the average V-band 
dust attenuation Ay. If we divide our sample into sources with 
dust attenuation lower than the sample median (i.e. Ay < 1.9) 
and sources with dust attenuation higher than the sample me¬ 
dian (Ay > 1.9), we hnd that galaxies with the lowest Ay have 
a median dust temperature Tjust = 38 ± 2 K and a median dust 
luminosity Ldust = (2.2 ±0.8) x lO'^ Lq, while the highest Ay 
galaxies have a median dust temperature Ldust = 42 ± 2 K and 
dust luminosity Ldust = (4.6 ±0.8) x 10'^ Lq. While this is a 
very tentative result given the large error bars on our Ay and 


Ldust estimates, it may indicate that, in more dust-obscured 
galaxies, dust is more effectively heated, presumably because 
of a more compact distribution of dust, i.e. the higher Ay, the 
more concentrated dust grains are around the (young) stars, 
and so they s ee a stronger radiation held (we return to this 
discussion in ^7.11) . A better sampling of the SEDs (possibly 
combined with known redshifts) is needed to reduce the sta¬ 
tistical error bars on these parameters in order to investigate 
this correlation in more detail. 

7. DISCUSSION 

7.1. The average spectral energy distribution of ALESS 
SMGs 

In this section we analyze the best-ht MAGPHYS ultraviolet- 
to-radio SEDs of our ALESS SMGs. In Eig. [TT] we plot 
the best-ht model SEDs of all 99 main ALESS SMGs in the 
rest-frame (calculated using our photometric redshift for each 
SMG), in observed hux units, highlighting the range of intrin¬ 
sic SEDs of these sources. We obtain the ‘average ALESS 
SED’ (plotted in black in Eig. fTTI) by computing the simple 
average Hux of all the best-ht SEDs at each rest-frame 
wavelength (weighting the average by infrared luminosity or 
sub-millimeter hux does not change the result). We also plot, 
for comparison, the SEDs of two prototypical local starburst 
galaxies, Arp220 (a ULIRG) and M82 (a d warf starburst), 
the ave rage SED of 16 local ULIRGs from Ida Cunha et 51 
(I2010al) . and the z = 2.3 sub- millimeter galaxy SM MJ2135- 
0102 (the ‘Cosmic Eyelash’: ISwinbank et al.ll2010l) . all nor¬ 
malized such that their total infrared luminosity is the same 
as the median total infrared luminosity of the whole sam¬ 
ple i.e. 3.5 X 10*2 Lq, at the median redshift z = 2.7. The 
spectral energy distribution of the average ALESS SMG in 
the infrared (A > 5^m) follows closely that of Arp220 (with 
the exception oT the 9.8-^m silicate absorption feature that 
is included in the Arp220 template but not in the MAGPHYS 
models). In the rest-frame ultraviolet to near-infrared range, 
both the Arp220 template and the average local ULIRG tem¬ 
plate are fainter and redder than the average ALESS SMG, 
which could be due to a different stellar content and/or dust 
attenuation. This is consistent with high-redshift SMGs being 
more massive and/or having a more extended dust distribu¬ 
tion than local ULIRGs. While Arp220 may provide a suit¬ 
able template for the (rest-frame infrared) dust emission of 
high-redshift SMGs, caution must be taken when using this 
template to interpret or extrapolate their (rest-frame UV-to- 
near-IR) stellar emission. Eig.[TT]also shows that the average 
emission by ALESS SMGs is very different from the starburst 
M82, which shows hotter dust emission, signihcantly higher 
optical/near-infrared Hux, and a redder stellar continuum. It is 
clear from this comparison that the full SED of high-redshift 
SMGs is not simply a scaled-up version of a lo cal starburst. 
Interestingly, iMenendez-Delmestre et al.l (120091) hnd that the 
mid-infrared spectral properties of SMGs are similar to those 
of local starbursts, but due to the poor sampling of the rest- 
frame mid-infrared emission, we cannot verify whether this is 
also the case for the ALESS SMGs. 

In Eig. [TT] we compare the average SEDs of our 77 
optically-bright sources and our 22 optically-faint sources. 
Not surprisingly, the ultraviolet-to-near-infrared emission of 
the optically-faint sample is signihcantly fainter (about an or¬ 
der of magnitude) than for the optically-bright sample, and the 
slope of the stellar continuum is also redder. This is consis¬ 
tent with the signihcantly higher inferred Ay of these SMGs 
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Figure 11. Best-fit model spectral energy distributions of the 99 ALESS sources obtained from our fits. In grey, we plot the SED of each individual source (we 
use our photometric redshift estimate to plot each SED in the rest-frame). We also plot the average SED obtained by averaging the flux Fu of each of the 99 
best-fit mode l SEDs at each (re st-frame) wavelength. For comparison, we pl ot the template SEDs of the ‘prototypical’ local ULIRG Aip220, the local starburst 
galaxy M82 iSilva et a T]|199j8i), the average SED of 16 local ULIRGs from Ida Cunha et all <2010al) . and the ^ = 2.3 sub-millimeter galaxy SMMJ2135-0102 
(‘Cosmic Eyelash’ dSwinbank et al.ll2Q10n . These four templates are normalized such that their total infrared luminosity is the same as the median total infrared 
luminosity of the whole sample i.e. 3.5 x 10^^ Lq, at the median redshift z = 2.7. 



A / fim [rest —frome] 

Figure 12. Comparison between the average SED of the 77 optically-bright 
sources and the 22 optically-faint sources. The shaded regions show the 16th 
to 84th percentile range for each subsample. 

as discussed in ^5.11 Interestingly, the average infrared SED 
of these SMGs peaks at shorter wavelength than the optically- 
bright average SED, which could indicate that these sources 
have on average higher dust temperatures. This tentative ev¬ 
idence for warmer dust in the optically-faint sources, com¬ 
bined with higher inferred dust attenuation Ay may indicate 
that these galaxies are more compact than their optically- 
bright counterparts, which would increase the dust column 
(hence higher Ay) and provid e more effective d ust heating by 
the stellar radiation field (e.g. lElbaz et al.ll201 ll) . 

To understand the variation of intrinsic SEDs of the 99 
ALESS SMGs, we plot, in Eig.[T3] the average SEDs in dif¬ 
ferent bins of (a) redshift, (b) observed ALMA 870-^m flux. 


(c) average E-band dust attenuation Ay, and (d) and total dust 
luminosity. We note that the scatter within each bin is large 
and in some cases larger than the difference between the bins, 
but nevertheless this is a useful way of understanding average 
variations of SED shapes as a function of different properties. 
A larger sample with more complete sampling of the SEDs for 
all sources is needed to perform a more quantitative statistical 
analysis. 

Fig- flBl al shows that the average intrinsic SEDs of the 
ALESS SMGs are different depending on their redshift: as 
redshift increases, we select intrinsically brighter sources with 
that peak at lower wavelengths in the (rest-frame) infrared i.e. 
have warmer dust. Eig. fOl^ b') shows that the average intrin¬ 
sic SEDs of SMGs of different observed 870-^m flux do not 
vary signihcantly: this does not mean that the SEDs do not 
vary within each flux bin, but rather that there is no system¬ 
atic variation with sub-millimeter flux, i.e. we are not select¬ 
ing a particular type of SED in different flux bins. If all we 
know about an SMG is its 870-/rm flux and redshift, then the 
average ALESS SED is a fair approximation for its intrinsic 
SED regardless of its sub-millimeter flux. Eig. fOl' c'). illus¬ 
trates how the average dust attenuation is a crucial parameter 
driving the variability of optical SED shapes in the ALESS 
sample; as expected, as Ay increases, the stellar continuum of 
the galaxies becomes systematically fainter and redder. We 
hnd that the infrared emission in different Ay bins remains 
almost constant, with only a small increase of the total in¬ 
frared luminosity and dust temperature towards higher Ay. 
The fact that the infrared SEDs vary little with Ay is a pos¬ 
sible indication that the main difference between galaxies in 
these bins is viewing angle (i.e. galaxies with higher Ay being 
observed more edge-on), since viewing angle does not affect 
the shape and normalization of the dust emission because it 
is isotropic (for an illustration of the effect of viewing angle 



























The physical properties of ALESS SMGs 


17 


redshift 



X / [rest-frome] 



X / jum [rest-frome] 




Figure 13. Variation of the intrinsic SED shapes of the 99 ALESS SMGs with different properties: (a) redshift; (b) observed 870-/rm flux measured with ALMA 
(in mJy); (c) average V-band dust attenuation; (d) total dust infrared luminosity (in Lq). In each panel, we plot the average SED in four different bins of the 
parameter being considered, indicated in the top right-hand com er (the number in brackets indicates the number of sources falling in each bin). For reference, we 
also plot the average SED of all the sources (shown in Fig. II H all as a black dotted line. All the SEDs are normalized to the average ALMA 870/xm flux of the 
sample, 4 mJy, at rest-frame 250^m, which corresponds to observed wavelength 870/rm for a source at z = 2.5 (indicated by the black circle). 


on SED s computed using a radiative transfer code, see fig¬ 
ure 8 of LTonsson et al.ll2010l) . The small observed variation 
in the shape and normalization of the dust emission in differ¬ 
ent bins (in the sense that higher Ay sources are on average 
slightly hotter and more infrared-luminous) could be due to 
variations on the compactness of the sources, which as dis¬ 
cussed above can be driving the difference between the SEDs 
of optically-bright and optically-faint sources. In this sce¬ 
nario, in more compact galaxies, the dust column is higher 
(hence higher Ay), and the dust gets more effectively heated 
because it is closer to the heating source (the stars) and hence 
feels a stronger radiation field. In Eig. fOl' d'). we plot the av¬ 
erage SEDs in bins of total dust luminosity. The peak of the 
dust emission shifts to lower wavelengths as the dust luminos¬ 
ity increases, whic h re fiects the luminosity-temperature rela¬ 
tion discussed in ^6.31 We also note that the average stellar 
continuum gets redder as the dust luminosity increases. This 
is mainly caused by the fact that the highest dust luminosity 
sources have higher stellar masses (and thus stronger near- 
infrared stellar bumps) and also, to a lesser extent, because 
they have slightly higher average Ay. 

7.2. The nature of SMGs: comparison with the ‘star-forming 
main sequence ’ 


In this section we discuss, based on our results for the 
physical properties of the ALESS SMGs, how they compare 
with the ‘normal’ galaxy population at their redshifts. The 
place occupied by SMGs in the star formation rate vs. stel¬ 
lar mass plane is a subject of great ongoing discussion in 
the literature (e.g.jHainline et all TOlU : iMagnelli et al1l2012l : 

iMichalowski et al.ll2012h . Some studies argue that SMGs are 

extreme starbursts forming stars at much higher rates than 
other galaxies of the same stellar mass at the same redshift, 
implying that these intense starbursts ha ve to be powered 
by extreme events such as major merger ([P addi et al l 120071 : 
lEngel et 51120101 : Hainline et alJl20flt Magnelli et al.ll20f^ . 

while some cosmological simulations establish SMGs simply 
as the high-mass end of the normal galaxy population, imply¬ 
ing that whatever ‘secular’ processes govern the star forma¬ 
tion histories of fainter gala xies are also at p lay in the sub¬ 
millimeter bright population (iDave et al.ll2010h . 

Our careful analysis of the star formation rate and stellar 
mass of the ALESS sample (a complete and unbiased SMG 
sample) gives us a unique opportunity to address this ques¬ 
tion. In order to compare the ALESS SMGs with the gen¬ 
eral galaxy population at their redshifts, in Eig. [14] we plot 
the star formation rates against stellar masses of our sources 
in two redshift bins (1.5 < z < 2.5 and 2.5 < z < 4.5), and 
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Figure 14. Comparison of the stellar masses and star formation rates of our ALESS SMGs with the ‘star-forming main sequence’ in two different redshift bins: 
1.5 < z < 2.5 (left) and 2.5 < z < 4.5 (right). Each galaxy is colour-coded by its 870-/rm flux density measured from the ALMA observations (Hodge_et^ 

S The dark grey solid lines show the position of the star-forming main sequence at each redshift z = 2 (left) and z = 3.5 (right) as given bv ISneagle et a . 

, based on a compilation of different studies in the literature spanning a wide range in stellar masses and redshifts. The star formation rate of main sequence 
galaxies is given by: log(SERMs(z, 4 f*)/M 0 yr“*) = (0.84 —0.026f(z))log(M* /Mq)-f0.1 lt(z) —6.51, where /(z) is the age of the Universe at redshift z i n Gyn . 
The grey dotted lines indicate a factor of 3 above and below this main sequence. Eor reference, the black star shows the ‘average SMG’ of IDaddi et al.l 1200711 
(with log(M* /Mq) = 10.6 and log(SER/M 0 yr“*) = 3). The black open circle in each panel shows the median stellar mass and star formation rate of the ALESS 
sources in that redshift range (with the error bar indicating the standard deviation). 


compare them with the ‘star-forming main sequence’ (i.e. the 
observed correlation between the stellar mass and the star for¬ 
mation rate of the general galaxy population) at those red- 
shifts. The exact slope, normalisation and redshift evolution 
of the main sequence still have large uncertainties and vary 
be tween studies; here w e chose the main sequence definition 
of ISneagle et al.l (120141) . which is based on a compilation of 
different studies in the literature spanning a wide range of stel¬ 
lar masses and redshifts. We plot their main sequence relation 
at z = 2.0 and z = 3.5 as a grey solid line in the left- and right- 
hand panels of Fig. [14] respectively. 

Fig.[T4lshows that, contrary to less-active galaxies in mass- 
selected samples, there is no strong correlation between the 
stellar masses and star formation rates of the AFESS sources 
at fixed redshift; this is because the sub-millimeter selection 
essentially selects in star formation rate. We find that, at 
z ~ 2, our SMGs have on average M* ~ 4.5 x lO'® Mq and 
SFR~ 205 Mq yr“*, which implies that their specific star for¬ 
mation rate is about 0.7 de x lowe r than that of the ‘average 
z = 2 SMG’ of iDaddi et al.l (120071) (plotted for comparison). 
With the relatively more modest specific SFRs of our SMGs, 
we find that 17 zb 2 (i.e.^49%) of them lie significantly above 
the main sequence at that redshift (i.e. they have star for¬ 
mation rates over a factor of 3 higher than ‘main sequence’ 
galaxies of the same stellar mass), while 18 zb 2 (~ 51%) are 
consistent with being in the high-mass end of the main se¬ 
quence. At higher redshift (2.5 < z < 4.5), the star formation 
rates and stellar masses of the SMGs are typically higher than 
at z — 2, with M* ~ 10" Mq and SFR~ 540 Mq yr“' on av¬ 
erage. However, since the speci fic star formation rate of main 
sequence galaxies predicted by ISneagle et al.l (12014l) contin¬ 
ues increasing with redshift, this means that a smaller frac¬ 


tion of our SMGs lie significantly above the main sequence at 
those redshifts. We find that only 11 zb 2 SMGs (i.e.~ 27%) 
have SFRs at least 3 times higher than main sequence galaxies 
of the same stellar mass at the same redshift, with most SMGs 
(73%) being consistent with the main sequence. Our results 
depend on how well the evolution of the star-forming main se¬ 
quence with redshift is constrained. If the specific star forma¬ 
tion of main sequence galaxies plateaus at z > 2 as suggested 
by so me studies (e.g. iWeinmann et al.ll201 ll; iGonzalez et al.l 
|20I3), then the fraction of high-redshift SMGs that are sig¬ 
nificant outliers would increase. The wide stellar masses and 
positions with respect to the star-forming main sequence of 
the AFESS SMGs may be an indication that th ese galaxies are 
not a uniform population (as suggested by e.g. lHavward et all 

l2ml . 

8 . SUMMARY & CONCLUSIONS 

In this paper we have presented an update of the MAGPHYS 
code that allows us to fit simultaneously the ultraviolet-to- 
radio SEDs of an unbiased sample of SMGs from the AFESS 
survey. This new version of the code allows us to constrain 
photometric redshifts simultaneously with all other physical 
properties, and explores a realistic parameter space for high- 
redshift star-forming galaxies, including complex star forma¬ 
tion histories, and a wide range of dust attenuation and emis¬ 
sion parameters to include different possible dust contents and 
distributions. Our main results are summarized as follows. 

• Our redshifts are consistent with the classic method of 
deriving photometric redshifts by fitting only the stel¬ 
lar emission (i.e. no dust emission). Our method im¬ 
proves the redshift constraints in the case of optically- 
faint sources for which optical observations are limited. 
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by including information from infrared and radio wave¬ 
lengths. An advantage of our method is that the un¬ 
certainty in photometric redshift is naturally included 
in the uncertainty in the intrinsic physical parameters 
of the galaxies such as the stellar mass, dust attenua¬ 
tion, dust luminosity, dust temperature and star forma¬ 
tion rate (and vice-versa). 

• We derive the median properties of the full sample of 99 
ALESS SMGs by has a median redshift z = 2.7 ±0.1, 
has a median stellar mass of ~ 8.9 ±0.1 x 10'** Mq, a 
star formation rate of ~ 300 Mq yr“', average V-band 
dust attenuation Ay ~ 2, total dust luminosity Ldust =— 
4 X 10'^ Lq, total dust mass Mdust — 6 x 10^ Mq, 
and luminosity-averaged dust temperature Tdust — 43 K. 
These physical properti es are very similar to t he prop¬ 
erties of local ULIRGs (Ida Cunha et aLll^lOah . 

• We find that the 22 optically-faint ALESS SMGs are 
likely to be at higher redshifts than the optically-bright 
sources, with a median photometric redshift of Zphot = 
3.7. Our analysis of the likelihood distributions of the 
physical parameters of these galaxies indicates that they 
are consistent with having similar overall properties as 
the bright sources, except for significantly higher dust 
attenuation. This could indicate that these SMGs are a 
population of high-redshift SMGs that are either very 
compact or very edge-on compared to sources of simi¬ 
lar sub-millimeter fluxes that are more optically-bright. 

• Using our multi-wavelength SED fits, we provide new 
SMG templates that should be more appropriate to in¬ 
terpret and/or extrapolate SMG observations than cur¬ 
rent local galaxy templates. 

• When we compare the star formation rates and stel¬ 
lar masses of the ALESS SMGs with the star forma¬ 
tion main sequence at that redshift, we hnd that, at 
z ~ 2, about half of the SMGs lie significantly above 
the main sequence (with SERs 3 times higher than main 
sequence galaxies of the same stellar mass), while the 
other half falls in the high-mass end of the main se¬ 
quence. At higher redshifts (z — 3.5), the SMGs tend 
to have higher SERs and stellar masses than at z ~ 2, 
but if we include the evo lution of the star-form ing main 
sequence with redshift (ISneagle et al.ll20f^ , we find 
that only about a third of the SMGs can be consid¬ 
ered extreme, high-SER outliers. This suggests that the 
ALESS SMGs are not a uniform population, including 
galaxies that may be extreme starbursts but also galax¬ 
ies with SERs that are similar to those of the general 
population of galaxies at the same redshift. However, 
in order to fully understand what this means in terms 
of the star formation mode of SMGs, we need to more 
detailed information such as spatially-resolved imaging 
and spectroscopy (to determine, for example, their dy¬ 
namics and merging state), and we also need to under¬ 
stand the processes that drive the star-forming main se¬ 
quence at high redshifts in more detail. 

Our method is designed to optimally extract as much in¬ 
formation as possible on the physical properties of SMGs 
(or, more generally, any sample of high-redshift star-forming 
galaxies with a wide range of dust contents and properties) us¬ 
ing integrated multi-wavelength observations. We will make 


the new MAGPHYS model library, as well as our average SMG 
templates, available to the community on the MAGPHYS web¬ 
site (www.iap.fr/magphys). 
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Table A1 

Median-likelihood estimates (and confidence ranges) of several physical parameters for each ALESS SMG 


ID 

2phot 

log(M./MQ) 

log(SFR/M0yr*) 

log(age*j/yr) 

Av 

log(M*/LH) 

log(Ldust/C0) 

Tiu . t/K 

log(Mdust/M0) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(V) 

(8) 

(9) 

(10) 

ALESSOOI.I 

4.78!?:“ 

10.97^:^ 

2.78^:« 

8-29!“:S 

2.0+l-‘] 

-0.12^:^“ 

12.82!“:« 

42+20 

^^-14 

9 19+0.49 
^••^^-0.46 

ALESS00I.2 

5-22!“* 

io.9o^:‘3 

2.62^:« 

8-36!“:“8 

1-2^1 

-O.lStJ? 

1 2 66+*^-^^ 
1Z.OD_q5-7 

45+21 

^^-16 

8.69!“:“* 

ALESS001.3 

T 1-^ + 1.10 

■^•^^-1.85 


2 46+0-"^2 
^•^‘^-0.92 

8.26!“:’“ 

3.0!]:’ 

0 OS+0-61 
^•^^-0.52 

12 55+0-39 
iZ.jj_o 88 

44!f^ 

0 ':5q+0.44 

°-7^-0.24 

ALESS002.1 

3.i2^:« 

I0.72tg] 

3.08t““ 

7 ^5+0.81 
'•^°-0.14 

2 9+0.7 
^•^-0.8 

-0-53ti 

13.02!“:*^ 

56+^ 

JO_24 

8.68!“:““ 

ALESS002.2 


1 1 0O+0-A9 

il.UU_0 5i 

2-77^:l“ 

0 94+0.53 
*••^^-0.68 

9 5+1.4 
”^•^-1.5 

0.19!“:^ 

12.84!“:““ 

48!]* 

8.64!“:““ 

ALESS003.1 

a qo+0.90 

11 37+0- 

^^•■^'-0.51 

2 93+0.29 
^•^°-0.40 

8 cn+0.34 
-0.70 

7.U-09 

0 ^9+0.30 
-0.56 

13.06!“;““ 

43!‘“ 

8.99!“:“;; 

ALESS005.1 

3 57+0.05 
■^•^'-0.20 

11.46^:“ 

3 10+O-02 

Q 91+0.30 
-0.18 

^•+-0.1 

-0 50^- 

‘^•'^^-0.17 

1 9 99+0.06 

^^•7'^-0.05 

40!* 

9.08!“;““ 

ALESS006.1 

0 QQ+0.15 
^•^°-0.10 

io.83^:«] 

2 92+0.04 
^•^^-0.07 

7 34+0.44 

'•°^-o.oi 

^•^-0.2 

-0 73+^-3^ 
^•'^-0.02 

12 75+0-05 
' -0.03 

43!“ 

9 17+0.02 
^••^'-0.20 

ALESS007.1 

I.98!<J;“ 

11 36^06 
ii.OO_Q_24 

9 91 +0.45 
^•■^^-0.28 

^--^^-0.45 

^•^-0.7 

-0.02^:]“ 

12 59+^-^^ 
i^.j7-o 28 

35!** 

^•^0-0.19 

ALESS009.1 


ii.75^:«» 

3 15+0.27 
•^•^^-0.66 

8 55+0.39 
-0.37 

2-51,:? 

o.o6!“:« 

13-24!“:“‘ 

53+13 

J7-19 

0 05+0.40 

0-0^-013 

ALESSOIO.l 

^•^^-0.01 

10.08!];;!]] 

9 78+0.53 

^•^O-003 

7 71+0.53 
-0.00 

3.2!“:“ 

-0 13+0-03 
^•■‘•^-0.24 

12-28!“:i 

4I!“ 

9.09!“:“* 

ALESSOIl.l 

3.03!<;:« 

11 45+0.10 
i^-"tO_0 85 

2.73!“:““ 

8 79+0.26 
' -0.62 

9 9+0.7 
■^••^-0.8 

^•^^-0.47 

12-84!“:]1 

45!“, 

9 99+0.07 
^•^^-0.24 

ALESS0I3.1 

3-22U;:i 

10 Q2+0-11 
lU.yz_o_28 

2.62!“:““ 

8 77+0.44 
-0.39 

1 8+0-3 
^•°-0.7 

-0 16+^'^3 
'-'•^'^-0.30 

12 66+^-^^ 
iZ.OD_Q 27 

42!* 

9 39+0.19 
^••7^-0.22 

ALESS0I4.1 

a OO+0.24 
•^-^°-0.26 

JO 95+0-13 

3.08!“:*“ 

7.98!“:““ 

7 7+0.4 
^•'-0.9 

-0.22!“:““ 

13 06+^-^^ 

i-7-UD_Q 13 

44!i4 

8.98!“:“* 

ALESS0I5.1 


1 1 49+0.15 
^^•^^-0.39 

2 4«:+0.34 
^•^'^-0.59 

8.89!“:““ 

3.4!]:t 

0.38!«:« 

12 64+^-^^ 

iZ.04_o 33 

34!‘“ 

9 29+0.18 

^••7^-0.20 

ALESSOI5.3 

7} 49+4.11 
•^•^^-2.34 

10 43+0-1'7 
^^•^+-0.63 

7 11 +0.68 

^•^^-1.43 

3 43+0.79 
-0.78 

1.3!“:“ 

-0.13!“:““ 

12 15+0-56 
iZ.iJ_i 20 

42+19 

^^-16 

Q 51 +0.54 
-1.00 

ALESS0I7.1 

9 19+0.35 
^•^^-0.60 

11 37+0-01 
-0.55 

2 46+0-21 
^•^^-0.16 

3 79+0.20 
^ -' -0.27 

^•^-0.4 

0.04!“:“* 

iz.oU q 26 

32+4 

■74-0 

9 'to+0.14 

^••7-7-0.20 

ALESS0I8.1 

^.UJ_020 

11.89!°;“ 

2 40+0.38 
^•^^-0.14 

g 1 9+0.00 
^••^^-0.50 

^•^-0.7 

0 25+0-00 
^•^'^-0.30 

12 69+^-^^ 
iz.oy_o.09 

37!“ 

8.99!“:^ 

ALESS0I9.1 

9 CT+O.OO 
•^•'^•^-0.56 

11 23+0-21 
^^•^•^-0.66 

2 94+0.34 
-0.38 

8.44!“:« 

3.3!“:“ 

0 ^5+0.36 
-0.68 

13.01!“:“’ 

59+14 

JV_12 

8.68!“:““ 

ALESS0I9.2 

2.17^:36 

10 40+0-12 
^^•^^-0.49 

i-8i!!]:i^ 

8-68!“:^* 

i.o!“:f 

-0 16+^-^^ 
-0.28 

11.86!“:“’ 

34!*“ 

8.64!“:““ 

ALESS022.1 

9 49+0.46 
^•^^-0.75 

11.60!“:“ 

2.66!“:““ 

8.90!«:“f 

2-i!!]j 

0 1 1+0.22 
^•^^-0.18 

12 77+0-31 
' '-0.42 

40!f 

Q 05+0.25 
°-°^-0.23 

ALESS023.1 

4 Q7+1.55 
^•^^-0.95 

ii.i8!“:i“ 

2 9«:+0.36 
^•^^-0.45 

8 78+0.45 
°-^^-0.76 

2>.J-i 4 

0 15+0.43 
^•^'^-0.64 

13 00+^-^^ 

i7.UU_Q 34 

49!]’ 

8.90!“:““ 

ALESS023.7 

3.08;^;f^ 

10.28!“:i 

1 95+0.70 
^•^*^-1.08 

0 4C+O.72 

o-^J-0.81 

1.3!“:“ 

-0 16+^-^^ 
'-'•^‘^-0.52 

12 on+^-^2 

iZ.UU_Q 82 

42!]* 

8.47!“:“’ 

ALESS025.1 

2 57+f>'50 
^•'^'-1.25 

10.53!“:’“ 

2 04+0.32 
^•*^-0.50 

0 47+0.44 
*•^'-0.56 

2.2!“:“ 

-0.11!“:“^ 

12 96+^-^^ 

iZ.VD_Q 53 

48!^* 

8.96!“:“* 

ALESS029.1 

3 52+0-95 
•^•'^^-0.54 

11.48!“:*“ 

2 97+0.35 
^•^'-0.39 

8 57+0.35 
-0.70 

2 9+1.1 
^•^-1.0 

o.i4!“:“f 

13 06+^-^^ 

i-7-UD_Q 21 

49!1J 

8.76!“:]’ 

ALESS031.1 

4.22!] 

1 1 45+0.12 
^^•^^-0.51 

2 92+0.39 
^•^^-0.49 

8 5^+0.33 
"■^^-0.65 

7 8+1-0 

^•°-0.9 

0 1 1+0-33 
^•^^-0.54 

13 01+^-^^ 

i7.Ul_o 37 

44!j| 

8.98!“:“’ 

ALESS035.1 

9 co+0.95 
•^•'^°-0.86 

11 24+0-20 
^^•^^-0.56 

2 04+0.37 
^•^^-0.54 

8 51+0.46 
-0.82 

2.6!]:’ 

0 07+0-54 
-0.62 

1 2 92+0-30 

46!]* 

8.64!“:““ 

ALESS035.2 

+57!i:?? 

10 21+0-23 
^^•^•^-0.55 

7 17+0.51 
^•^^-0.61 

8.26!“:““ 

4.9!^:* 

0 42+0.87 
-0.82 

12 17+0-47 
^^••^'-0.48 

43+1^ 

^■7-15 

8.40!“:““ 

ALESS037.1 

2.72!«:i| 

11 30+0-02 
li-7U_0,30 

2-80!“:“] 

8.51!“:““ 

i-8!!]j 

-0 12+0-14 
^••^^-0.34 

12.84!“;““ 

53+14 

^■7-16 

0 49+0.25 
*•^^-0.15 

ALESS037.2 

a o't+0.64 
•^•°-^-0.55 

JO 55+0.10 

7 99+0.52 
^•■^•^-0.44 

8-27!“:S 

o.8!“:l 

-0 39+0-25 
'^-7^-0.30 

1 2 31+0-50 
^^••71-0.45 

47+19 

^'-16 

0 94+0.52 
*-7^-0.47 

ALESS039.1 

9 99+0.50 
^•■^•^-0.35 

10 55+0.09 

2 94+0.36 
-0.20 

8 77+0.44 
°-++-0.35 

1 4+0.3 
^ -^-0.4 

-0.28!“:** 

12 31+0-35 

iZ.71-0 13 

32!f 

9 99+0.13 

^•17-7-0.13 

ALESS041.1 

9 17+0.61 
^•^'-0.65 

11.35!“:““ 

2 25+0.62 
^•^^-0.49 

3 9J+0.38 
^•^^-0.68 

1 7+^-2 
'-1.2 

o.o8!«:g 

12 43+^-'^"^ 

iZ.H7_o_27 

34!f 

9 06+^-^^ 
y.UO-0 26 

ALESS041.3 

2.97!4;46 

9.86!“:“> 

9 19+0.84 
^•■‘•^-2.70 

8 95+1.13 
^••^'-’-1.08 

0 7+2-^ 
^•^-0.6 

—0 24+0-80 
^•^^-0.60 

12 11+0-71 
-11-2.46 

42+19 

^^-16 

8.42!“:’“ 

ALESS043.1 

9 no+0.50 
^•^°-0.60 

11 31+0-09 
ii.j*i_0 2o 

1 90+0-56 
^•^^-0.39 

9 07+0.28 
^•^'-0.46 

7 8+1-0 
^•°-0.8 

0.33!“:““ 

1 9 1 C+O.40 
14.17-0 21 

34!*“ 

Q 77+0.26 

'^-0.30 

ALESS045.1 

9 17+0.86 
•^•^'-0.59 

11 71+0.11 
^^•'^-0.26 

2 51+0.39 
^•^^-0.46 

8 8^+0-30 
°-°^-0.41 

2-8!^:? 

o-28!!j:i 

12.77!“:““ 

38!‘“ 

0 99+0.22 
^-7-0.22 

ALESS049.1 

7 07+0.14 
^•°-^-0.05 

10.58!“:*“ 

7 8'3+0.09 
^•°-^-0.05 

7 75+0.83 
/• /J-0.19 

1.9!“:“ 

-0 33+0-10 
^•■7-7-0.31 

12.83!“:“^ 

46+“ 

9 94+0.04 
^•^^-0.12 

ALESS049.2 

2 57+0.66 

^•^'-0.25 

11 02+0-15 
^^•^^-0.36 

2 C4+O.39 
-0.25 

8.60!“:“ 

^•^-0.6 

-0 11+0-18 
^••^^-0.42 

12 60+^-^^ 
iz.ou_o 19 

51!‘“ 

Q 70+0.21 
^•^^-0.21 

ALESS051.1 

1 99+0.19 
^•■^•^-0.10 

11 17+0-07 
-0.24 

1 Q7+0.20 

-0.33 

9-05!“:“’ 

^•^-0.6 

0 11+0-17 

^•^i-0.19 

12 13+^-^^ 

1Z.17_q 12 

31!‘“ 

9 ^^+0.09 

ALESS055.1 

7 70+0.25 
^•"^°-0.20 

9 97+0.27 
-0.21 

9 99+0.18 
-0.22 

8.00!“:“ 

1 1+0-8 
^•^-1.0 

-0 39+0-16 
^•-'^-0.33 

12 15+^-^^ 

29!‘“ 

9.05!“;]“ 

ALESS055.2 

4.68!2;S9 

10-42!“:““ 

2-36!“:^* 

8.24!“:““ 

5 4+3-7 
-3.4 

-0.80 

12 40+^-'^^ 

iz.hU-o_47 

44+18 

^^-16 

8.61!“:““ 

ALESS055.5 

9 99+0.14 
^•■^•^-1.06 

10-15!“:“ 

0.80!*:““ 

9-i3!?:^i 

0 0+^31 
^•^-0.0 

-o-i5!]]:§? 

10 19+1-53 
-0.81 

44!““ 

^■38!a38 

ALESS057.1 

1 9S+O.6O 
^•^°-0.50 

9 99+0.27 
^•^^-0.13 

9 99+0.41 
^•■^^-0.18 

7 70+^-*^^ 
'^-0.27 

2.4!“:“ 

-0.48!“:^“ 

1 7 '37+0.40 

IZ.J/-020 

42+15 

^4-14 

^•^^-0.42 

ALESS059.2 

i.48!<;:if 

9.76!“:*“ 

1 87+0.12 
-0.23 

7 97+0.76 

'•^'-0.30 

2 5+0.5 
Z.7-i_2 

-0 16+^'^^ 

-0 32 

11.88!“:*“ 

34!‘“ 

0 77+0.16 

° - '-7-0.22 

ALESS061.1 


10.58!“:“ 

9 14+0.14 
•^•^^-0.30 

7 97+0.89 
'•^'-0.00 

o-5!!J:i 

—0 29+*^'^*^ 
*^•"^^-0.71 

13.00!“;*“ 

C7+9 

J/-14 

0 57+0.62 
O.j'-o_22 

ALESS063.1 

2.08!<;;20 

11 20+001 
^^•^^-0.19 

1 93+0.00 
^•^°-0.34 

9 01 

^•^^-0.63 

1-2!!]:^ 

0 03+^-^^ 
^•^■7-0.14 

12.14!“:“* 

37+0 
•7'-13 

9 2f;+o.i3 
^•^'7_qo7 

ALESS065.1 

5.68!>:^9 

10 74+0-12 

iU./H_0 4g 

2 54+0-44 
^•^^-0.58 

8 79+0.47 
^•^■^-0.65 

1 9+1.8 
^•^-1.1 

-0.17!“:“^ 

iz.oo_Q_42 

44+17 

^^-16 

0 91+0.56 
^•^^-0.48 

ALESS066.1 

I.98!?:;;]] 

10.07!“:““ 

2-6l!!]j^ 

-0.34 

0 5^^'3 

^••^-0.3 

-0.80!“:^* 

12 44+0-34 
^^•^^-0.45 

49+10 

^^-10 

0 71+0.12 

6. /1-0 27 

ALESS067.1 

2.o8!<;:3o 

11.38!“:““ 

2.57!“:“] 

8.91!“:““ 

1 7+0.3 
^•'-0.7 

0 02+^-^^ 
^•^^-0.46 

12 70+^-^5 
iz./U-o_24 

38!“ 

0 07+0.19 
-0.13 

ALESS067.2 

1 C9+I.O6 
^•-^^-0.25 

Q 81 +0.087 
^•°^-0.25 

1 8S+0-35 
^•°^-0.40 

8.01!“:“^ 

2.0!“:“ 

—0 44+0-38 
^•^^-0.14 

1 1 05+0.36 
11.62-0 28 

42!]“ 

8-62!“:““ 

ALESS068.1 

9 78+1.90 

10 97+0-19 
^^•^^-0.57 

9 51+0-44 
^•^^^-0.53 

8-48!“:“^ 

3.2!]:^ 

^•^'7_o64 

12 69+^-^^ 

14.0^-0 40 

47!]’ 

8.67!“:““ 
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Table A1 

(continued) 


ID 

2phot 

log(M./M0) 

log(SFR/M0yr') 

log(age^/yr) 

Ay 

log(M*/LH) 

log(Ldus,/Z,0) 

Tdust/^ 

log(Mdust/M0) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

ALESS069.1 

9 oo+O-SO 

1L43^:‘> 

9 9Q+0.58 
+•^^-0.70 

8 87+0-36 
°-°'-0.45 

2-7!i:J 

0 26+11-^^ 

12 50+0-4'7 

iz.::)u_Q5Q 

36+23 

2>0-io 

8.99“:“ 

ALESS069.2 


io-38U:i 

2.32“;5; 

0 99+0.59 
°-^-^-0.72 

4.2“:; 

A 90+0.96 

b-^o_078 

1 9 97+0.46 
^^•'^'-0.47 

43+18 

^■^-15 

8.64“:“ 

ALESS069.3 

4.6Stlil 

10 36^-^^ 

2.29“;i“ 

8 9/1+0.60 
°-^^-0.73 

^•'^-3.4 

0 49+0-81 
^•^^-0.80 

12 34+0-45 
iZ.9H-_o_48 

44+18 

^^-16 

Q cc+0.56 

ALESS070.1 

1 ^8+0.64 
i.JO_0 2o 

10.63!<>;“ 

2 65+0.14 

z.UJ-o 17 

8.01“;“ 

2.5“:° 

-0 47-^^-40 
^•^'-0.00 

12 63+0-21 

46“ 

8-84“:“ 

ALESS071.1 

1 77+0-21 
^-0.39 

1 1 30+0.18 
il.J*y_0 35 

9 7C+O.I4 
'-^-0.30 

8-80“;“ 

9 c+0.4 
^•^-0.5 

0 11+0-1”^ 
^•^^-0.37 

12.82““ 

43“ 

8-7i!;:i; 

ALESS071.3 

9 90+0.39 

^•^o_0 80 

Q ^4+0.132 
^•‘^^-0.42 

L62“;« 

8 1 8+0-78 
°-l°-0.83 

0.8“:; 

-0 37+11-30 
'-^••^'-0.35 

il6i“:;5 

37!“ 

8.27“:“ 

ALESS072.1 

c 09+1.65 
j.oZ_2 ig 

i^-^^-0.41 

2-24“;^l 

8 91 +0.43 
°-2l-0.58 


-0 2+0-^^ 
^•^-0.43 

12 77+0.38 
'-0.48 

44!" 

g 9:1+0.61 

o-^'7_q44 

ALESS073.1 

4 78+0.40 
^•'°-0.50 

10.64!;;;“ 

2 90+0-30 
^•^^-0.36 

-0.31 

-0.7 

-0 41+0-19 
^•^^-0.26 

12 89+ii-^^ 

iZ.O!^_035 

46+12 

‘tO-14 

8-97!^:;; 

ALESS074.1 

^•^^-0.24 

10.68!;;;2» 

9 17+0.45 

z,.i /-o_22 

8 87+0.45 
-1.07 

^•^-1.1 

o-i23:?4 

12.17““ 

32!" 

9.04“:® 

ALESS075.1 

1 08+0.24 

^•^°-o.io 

10.48!;;;“ 

2.65“;|3 

7 74+0.56 

^•'^-0.19 

^•'^-0.5 

-0-683:^; 

12.58“" 

39“ 

8-76“:“ 

ALESS075.4 

2.28:^'“ 

^•^o_0_73 

1.57“:^ 

8-0i!g;M 

0 5+111 
-0.4 

-0 44+0.31 
l^-‘l^-0.40 

11 52+0-31 

ii.9Z_i_27 

38!“ 

8-111:?; 

ALESS076.1 

21 07+1.71 
■^•^'-0.94 

11 ni+0.20 

^^•^^-0.49 

9 04+0.37 

^•°^-0.46 

8-3i!g:^I 

3 5+1.4 

2>.J-1_7 

A 17+0.47 

^•1 '-0.53 

12.91““ 

44!" 

8.90“:® 

ALESS079.1 

9 ^9+1.09 
■^•'^•^-0.86 

il.DO_Q3Q 

9 co+0.44 
^•^^-0.54 

8.76“;“ 

9 7+0.8 
^•^-0.9 

0 20+0-50 

12.78““ 

43!" 

8-66“:“ 

ALESS079.2 

1 68+0-01 
^•^°-o.oo 
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Note. — (1) Catalog ID; (2) photometric redshift; (3) stellar mass; (4) star formation rate; (5) mass-weighted age; (6) average V-band dust attenuation; (7) 
//-band mass-to-light ratio; ( 8 ) total dust luminosity; (9) luminosity-averaged dust temperature; (10) total dust mass. 







